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SECTION 1 


INTRODUCTION 


This report presents the derivation, description, and operating instructions for a 
computer program which measures the economic value of advanced technology 
features applied to long range commercial passenger aircraft. The program consists 
of three modules; an airplane sizing routine, a direct operating cost routine, and an 
airline return-on-investment routine. These modules are linked (as shown in Figure 
1-1) such that they may be operated sequentially or individually, with one routine 
generating the input for the next or with the option of externally specifying the input 
for either of the economic routines. 


The program, labeled TEKVAL, was written by General Dynamics/C on vair under 
NASA/LRC Contract No. NAS1-11343. The concept for TEKVAL grew out of Convair's 
participation in the Advanced Transport Technology (ATT) studies, also sponsored by 
NASA/LRC, in 1971-72. That study examined the potential benefits to 
commercial transports of supercritical aerodynamics, advanced propulsion with 
reduced noise, composite structure, and active control systems combined with 
reduced static stability. Configurations were optimized for maximum return-on- 
investment to an airline operating a fleet of such airplanes over a selected route 
structure; an intermediate step required the computation of direct operating cost as 
a function of distance for each configuration. 


During the course of the ATT study, a very simple airplane sizing technique was 
developed, based on the Brequet range equation. For this contract, that sizing 
technique has been greatly expanded and combined with the formerly separate DOC 
and ROI programs to produce TEKVAL. The derivation of TEKVAL draws heavily 
on ATT study results to derive factors which reflect the impact of the above mentioned 
advanced technologies as an integral feature of this program. 

The computer output of the SIZE routine consists of two pages. The first page lists 
the principal input specifications followed by derived dimensions and geometry and 
aerodynamic and performance parameters. The second page is a weight statement. 

The DOC module also provides two pages of output. The first is a detailed breakout of 
labor and material manufacturing costs corresponding to the above weights statement. 
The second page is a detailed listing of the various elements of the ATA direct opera- 
ting costs formula computed as a function of distance. In the DOC program, the user 
has the option of an alternate summary cost equation referenced to AMPR weight; 
when utilizing this option, the first output page is replaced by a summary cost state- 
ment at the head of the DOC output. 
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Figure 1-1. TEKVAl Program Flow 











The first page of output from the ROI section repeats a number of input constants, a 
few basic airplane design parameters, and a listing of block time, and includes DOC 
and passenger fares as functions of distance. The next page(s) shows a breakout of 
passenger traffic, load factor, operating costs, income, profit, and ROI for each 
city -pair of the input set of operating routes. This breakout is followed by a total 
system summary cf the same parameters. At the option of the user, the individual 
city pair listings may be suppressed. 

It should be noted that program input and output utilize the English system of units. 
However, for the purposes of this report all dimensions, including those associated 
with input and output parameters will be shown in terms of the scientific international 
system of units with the English equivalent in parenthesis. 


c 
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SECTION 2 

ANALYTIC DERIVATION 


2.1 AIRPLANE SIZING TECHNIQUE 

The airplane sizing technique is based on an inverse solution of the Brequet range 
equation. Airplane takeoff gross weight is derived as a function of payload weight and 
cruise conditions. Airplanes are sized for that corner of the range-payload curve 
(Figure 2-1) where maximum fuel capacity is required to carry the full passenger pay- 
load over the design range, constrained by a maximum takeoff weight limit. The 
assumed flight profile is shown in Figure 2-2. The airplane cruises at constant Mach 
number, continually increasing altitude to remain at the optimum lift-to-drag (L/D) 
ratio. Engines are sized to provide the necessary cruise thrust-to -weight ratio (T/W 
= D/L) at the initial cruise point. Takeoff performance is not computed; however, 
static thrust-to-weight provides a useful indication of expected field lengths via refer- 
ence to existing commercial aircraft performance. 


The Brequet range equation is written as 

R = In (W /W ) 

c' o 1 


(1) 


with 


R = Range km (nm) 

V^. = True airspeed km/hr (kt) 
L/D = Lift-to-drag ratio 


c» = Thrust specific fuel consumption kg/N-hr (lbm/lb-hr) 


= Initial cruise gross weight kg (lb) 
= Final cruise gross weight kg (lb). 


This equation, derived in most basic aerodynamics texts, is applicable only to cruise 
flight, the difference between W 0 and being the weight of fuel used in flying the dis- 
tance, R. In order to compute gross takeoff weight (GW), equation (1) is solved for 
the initial cruise weight, W 0 , to which is added takeoff and climb fuel weights (Wf fcQ 

and Wf , , respectively); 
cl 


GW = W 


exp 


Re' 

V t (L/D) 


+ W + W . 
to cl 


( 2 ) 
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Figure 2-1. Range - Payload Design Point 
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The range, R, is the design range less climb and descent distances. The final cruise 
weight, W~l , is the sum of operating weight empty, payload, and reserve and descent 
fuel weights. ' 

The program employs an empirical expression to generate an initial estimate of takeoff 
gross weight. A fuselage layout is derived from payload capacity and seating dimen- 
sions specified to the program. An airplane empty weight is computed from a large set 
of empirical relationships, using engines scaled to a given static thrust/weight ratio. 
Climb and descent performance are computed to give these distance and fuel increments 
Finally, a takeoff gross weight is extracted via equation (2). A single iteration is em- 
ployed for improved wing and engine sizing, and subsystem weight computations. The 
computation procedure is illustrated in Figure 2-3, which also constitutes a computer 
flow chart. The numbered paragraphs which follow correspond to the like numbered 
boxes in Figure 2-3. In general, simple geometric computations have not been listed 
in this derivation. Reference to the computer program listing and dictionary will clari- 
fy these operations. 

1. Fuselage Layout — The sizing process begins by generating a fuselage layout to 
enclose the required seating arrangement (see Figure 2-4). Floor dimensions of 
the coach section are determined first to set fuselage diameter; external diameter 
is made 10 percent greater than floor width to allow for wall thickness and lost 
space. First class section length is determined and one diameter is added to the 
length of these two sections to account for galley, wardrobe and lavatory space. 

The sum of these three sections make up the straight portion of the fuselage. The 
cockpit and tailcone areas are referred to as tapered sections; their combined 
length is estimated as a specified number of fuselage diameters. For area ruled 
fuselages, the length of the straight section is increased by 25 percent and an ex- 
tra 1/3 diameter is added to the length of the tapered sections: both of these con- 
stants were derived by reference to ATT study configuration layouts. The several 
equations used in this section are merely arithmetic manipulations of geometric 
input data concerning seating dimensions and arrangements. They can easily be 
followed by reference to the program listing and definition of terms. The princi- 
pal output values of this section are: 

D = Fuselage external diameter m (ft) 

L = Length of constant diameter section m (ft) 
st 

L = Length of tapered nose and tail sections m (ft) . 
tpr 

2. Initial Gross Weight Estimate — Empirical expressions were derived relating 
gross weight/payload weight ratio to design range (R) and Mach number for several 
potential technology combinations. The basic equation for gross weight / payload 
weight ratio is: 
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Figure 2-3. Sequence of Operations 
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LENGTH = LENGTH = NO. SEAT ROWS X 

1 DIAMETER SEAT PITCH 



DIAMETER = (NO. SEATS ABREAST X SEAT WIDTH + 
NO. AISLES X AISLE WIDTH) X 1. 10 


Figure 2-4. Fuselage Layout 



( 3 ) 


GW/W = (3 + 1.3 R/lOOO)k .. 
pi tech 

For design Mach numbers (M) greater than 0. 9, equation (3) is multiplied by the 
factor [l + 20 (M - 0.9) 2 ]. Aerodynamic and structural technology is reflected ' 
in the factor kteeh which has the following set of values: 

Aerodynamics Conventional Supercritical 

Structure 

Light Alloy 1.1 1.0 

Composite 1.0 0.9 


3 . Initial Engine Sizing — Representative baseline engine performance must be pro- 
vided to the model in the form of tables of net thrust and specific fuel consumption 
vs. Mach number and altitude. Also required are the engine sea level static 
thrust rating, basic weight,and nacelle diameter. All of these values pertain to 
an engine of a particular size. The program computes an engine scale factor by 
which all engine performance, weight,and dimensions are scaled up or down to 
match the thrust required at the design condition. 


An initial estimate of engine scale (ES) is made so as to give a static thrust-to- 
weight ratio as specified by input. Engine weight and nacelle diameter are scaled 
via the relations 


D/D = (ES) 0 * 5 
o 


(4) 


and 


1.11 

W/W = (ES) 
o 


(5) 


where D Q and W Q are the diameter and weight values for the 1. 0 scale engine. 
These relationships should hold over an engine scale range of 0.5 - 2. 0; beyond 
this range a different baseline engine should be selected. 


4. Wing Sizing — Wing area is computed from gross weight and wing loading. Wing 
dimensions of span, root and tip chords, and mean aerodynamic chord (m. a. c. ) are 
computed from area and specified aspect ratio and taper ratio (see program listing). 
Wing sweep angle is given at the 1/4 m. a. c. and is translated to the leading edge 
for use in computing wing thickness ratios. For conventional and supercritical 
airfoils, thickness/chord ratio (t/c) is related to design Mach number and leading 
edge sweep angle (Aj e ) by the expressions 
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conventional: t/c = 0.802 - M (Cos A, ) 

le 


0. 6775 


and 


supercritical: t/c = 0.896 - M(CosA^) 


0. 6775 


(6) 

(7) 


5. Tail Sizing — The longitudinal locations of wing and tail surface m. a. c. 's are 
given as fractions of total fuselage length. From these, tail arm lengths are 
determined. The tail volume coefficients of Figure 2-5 are then employed to 
compute horizontal and vertical tail exposed areas. In the case of a fin-mounted 
engine such as on the McDonnell -Douglas DC-10, vertical tail area is reduced by 
the projected area of the nacelle. Because of its relative length, the diameter of 
a fin-mounted nacelle is approximately 20 percent greater than a podded wing — 
or fuselage-mounted installation. This factor and a nacelle length/diameter ratio 
are combined to compute nacelle projected. 


6. Initial Cruise Lift/Drag Ratio (L/D) — This parameter is required for subsequent 
engine sizing; it also represents a good measure of the aerodynamic efficiency of 
the total airplane and is included in the output listing. Cruise L/D is computed as ■ 
the mathematically equivalent ratio of lift/drag coefficients (C^/C d ). 


Lift coefficient (C^ ) at initial cruise weight is first estimated by the equation 
C z = 0.95 (W/S)/ 1481 6 M 2 (8) 

where 

2 2 

W/S = Design takeoff wing loading N/m (lb/ft ) 
and 5 = Atmospheric pressure ratio at initial cruise altitude . 


The factor 0.95 is inserted as an estimate of the ratio of gross weights at initial 
cruise and takeoff conditions. For the second program iteration, this factor is 
modified by the results of the first pass. A table of pressure ratio vs. altitude 
yields the value of 6 for the inputted initial cruise altitude. 

The drag coefficient (C^) at initial cruise conditions is measured by the expression 
C d ■ C fe < A/S w> , < 9) 

with 

C = Equivalent skin friction coefficient (input) 
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Figure 2-5. Tail Volume Coefficients. 
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2 2 

A = Total airplane wetted area m (ft ) 

2 2 

and S = Wing area m (ft ) 
w ' ' 

Total airplane wetted area is the sum of the surface areas of the fuselage, wing, 
tail, and nacelles and pylons. Fuselage gross wetted area (Af^) is approximated 

by 


■= TTD (L + 0.7 L ) 
fus fus st tpr 


( 10 ) 


where fuselage dimensional parameters were specified early. Wing root fairings, 
or gloves, cover an appreciable portion of the fuselage side surfaces; these masked 
areas must be subtracted from gross area to prevent an overestimate of airplane 
drag. The cross-sectional area of an airfoil shape (A cs ) can be estimated by 

A = 0.7 c 2 (t/c) 
cs • 


where c is the airfoil chord and t its thickness. Figure 2-6 illustrates the geo- 
metric terms employed in defining wing glove dimensions. The wing root is 
assumed to intersect the fuselage in a plane canted 20 degrees as shown. The 
lateral distance from the airplane centerline to this intersection is labelled y ^ 
and is expressed as 


(D £us /2,COS(20 ’ ) 


= 0.47 D 


fus 


The ratio of the root chord of the glove (u in Figure 2-6) to the theoretical wing 
chord computed at (v in Figure 2-6) is termed Cg £ , and is inputted to the 
program. With the wing chord, v, computed from available geometric param- 
eters. Net fuselage wetted area can now be obtained as the gross area, from equa- 
tion (10), less the glove cross sectional areas: 

A = (2)(0.7) v 2 (t/c)c „ /Cos (20°) 
cs gJt 

= 1.49 v 2 (t/c) c ^ . (11) 

Cross sectional areas of the tail surfaces are neglected. 

For wetted area computations, the wing is divided into inboard and outboard sec- 
tions, with the break coming where the glove intersects the straight leading and 
trailing edges, as seen in Figure 2-6. The lateral location of the break point is 
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input as a fraction of the wing semispan. For the outboard section, wetted area 
(A 0 jj) is figured as 



c 

sw 



where 


( 12 ) 


c = Ratio of wing wetted area/planform area (input) 

sw 2 2 
S , = Planform area of outboard wing sections m (ft ). 
ob 

A value of 2. 02 is recommended for c sw . For the inboard section, equation (12) 
is modified by an expression involving Cg^: 



with 


c s -k 

sw lb 


0.043 (c , 
gl 



(13) 



Wetted area of inboard wing sections including gloves m (ft ) 


Theoretical planform area of inboard wing sections from y 
in Figure 2-6 m 2 (ft 2 ). 


gi 



Tail surface wetted areas are computed via equation (12), using the same input 
value of c 

sw 


Nacelle length/diameter ratio, (L/D)nac (input), and previously computed diam- 
eter are combined to compute a cylindrical nacelle wetted area. A 20 per- 
cent increment is added to account for pylons, giving the complete expression 


A = 1.2ttd 2 (L/D) 
nac nac nac 


(14) 


Equation (14) applies only to podded engine installations. For a fin-mounted en- 
gine, with its relatively long nacelle and increased diameter, 


A 

nac 


1.44 TiD 


nac 


(L/D). 


nac -fin 


(15) 


A separate (L/D ) ^ ^ is input for this configuration. 
Total wetted area is the sum of the above components. 
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7. Engine Resizing — Having computed L/D at initial cruise conditions, engine size 
is rescaled to give T/W = D/L. Required cruise thrust per engine is computed at 
weight W Q . The cruise thrust for the full-scale engine at initial cruise Mach and 
altitude is extracted from the appropriate table. Engine scale is recomputed as 
the ratio of these two thrust values. Engine weight and nacelle diameter are ad- 
justed to the new scale factor, as are wetted area and L/D. 

8. Fuel Volume —An initial estimate of total fuel quantity is required to determine 
fuel subsystem weights. In the Brequet formula, equation (1), the difference be- 
tween l\^and W^ is the weight of cruise fuel Wf cr . Range factor (RF) is the label 
applied to the ratio Vj. (L/D)/c ' . Equation (1) can be manipulated to produce the 
expression 

W fc/ W o = 1 " ex P (" R/RF). (16) 

An empirically derived increment of 0. 1 is added to this expression to account 
for climb, descent and reserve fuel. Total fuel weight (Wf tot ) is then obtained 
from 

Wftot = 0.95 GW [l.l - exp (- R/RF)] (17) 

where GW is still the current estimate of takeoff gross weight. 

Dividing (17) by fuel density then yields fuel total quantity in gallons. 

Total design range is inserted for R in evaluating equation (16). The range factor, 
RF, is computed at initial cruise conditions; L/D is already available, and c' is 
obtained via a table look-up. True airspeed (Vj.) is computed from the relation 

V = V /'fa 
t e 

= M(V/M)/Va' (18) 

where 

V g = Equivalent airspeed km/hr (kt) 
and a = Atmospheric density ratio. 

This form was chosen for convenience for all Vj. vs. M computations to minimize 
data storage requirements (it is used repeatedly in the climb and descent segments); 
the parameters V e /M and Va are both obtained from tables as functions of altitude. 

This entire fuel volume estimation step is omitted in the iterative pass. 
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9. Weight Empty — A series of empirical weight estimating relationships for struc- 
ture, propuls ion, and subsystems and equipment were derived from data on a 
large number of existing commercial transports. These equations are listed in 
Table 2-1. 

The impact of several advanced technology features (area ruling, composite struc- 
ture, active control systems, and supercritical aerodynamics) is accounted for by 
modifying factors applied to the basic equations. A series of three control factors 
are input to define which advanced technologies are to be employed. At this point, 
only the body (W^) and wing (W w ) structural weights are affected: 

Body: W b = W b (1 - 0. 016 Kar , (1-0.2 ^ MATL) . Q1 K ) (19) 

ACo 

where 


K = Fuselage area ruling factor (0 = without, 1 •= with) 

AK 

KtyLATL = < “' onstructi<m material factor (0 = light alloy, 1 = composites) 

K = Active control system factor (0 = without, 1 = with). 

ACS 

The above three factors are input constants. 


Wing: W = W • R . R • R A „ CT 

w w MACH MATL ACS 


( 20 ) 


For M > 0.85, R„ A = 1 + 0. 7 (M - 0.85). 

MACH 

FOT K MATL = R MATL = °- 885 - °- 00255 < b/CosA o/2 ,/t r- 
ForK AC S ' la " dK MATL = 0 - R AC S =0 ' 871 - 

te Vf 1 - VTl S1 ' Vs' 0 ' 912 ' 

The three R-factors are internally computed. If M < 0. 85 or or *"ACS 

are zero, the corresponding R-factors are set equal to 1. 0. 

The summation of all items listed in Table 2-1 yields weight empty (WE). 


10. Operating Weight Empty — The addition of useful load to WE gives OWE. Useful 
load weights consists of four items computed from the following equations; 
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Table 2-1. Weight Estimating Relationships 
(Output in Pounds) 


Structure 


Body: 


Wing: 


Tail: 


, 0. 3 0. 9 1 O 1 ^ 

W b = (0.2712 • GW • ULF) L f D 1 * Ut> 


W = 0.306 
w 


GW * ULF * S w • ^b/Cos A c/ / 2 )^ 


1000 (t/c) C 


0.62 


root 
0.9 


W, = 12. 8 • S, " (S. = exposed area) 

h, v h, v h, v 


Nacelle: W (ea) = 3846 * ES - 977 
nac 


Ldg Gear: W = 0. 046 • GW 

Xg 


Propulsion 


, 1.11 


Engine: W (ea) = W • ES 

eng eng 1 

Engine/Nacelle Sound Proofing (if employed): 


W (ea) = 0. 035 • RATING - 460 
sp 


Starting: W = 0.35 • W 

st eng 


0.65 


Controls: W = 120 
ec 


' [( L l + b/C0sA C / 2 )N en g /l “° 


0.294 


Water Injection Provisions (if employed): 

W . (ea) = 3.41 x 10 _3 • RATING 
wi 


Fuel System: 

Pumps; W = 1.1 X 10~ 3 • RATING /1.75 N + 0.266 N 2 

p l eng eng 


Distribution: W.. = 0.24 • N V RATING + 0.62Q 

dist eng 


0.7 
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Table 2-1 (Cont’d) 


Venting: 

W = 3.48 x lo" 

V 

’ 3 • N • RATING 

eng 


Control: 

W. = 1.116 • Q 0,5 

IC 


Refueling: 

W . = 4.9 • Q°* 
rt 

33 


Dump: 

W. = 0.159 • 

dmp 

q °.65 


0 75 

Sealant: W = 0. 282 • Q 

S X* 

Systems and Equipment 


Surface Controls: 

W = 340 + 
sc 

0.95 (S + S^ + S ) 
w h v 

-3 


Instruments: 

W. = 2.68 X 10 • GW + 165 

inst 


Hydraulics: 

W. = 0.48 
hyd 

(V 

0.45 (S + S u + S ) 
w h v 

. A V 1.06121 

b/C0E A c/z) 

1.3125 

i • + 

0. 849 

Avionics: 

W . = 1.3 • 
avi 

WAVI 


Electrical: 

W , =62.8 

/ew. . - W + W . \ 0 * 473 

elec 

+ 2 

\ fuel sys 6* 

• PAX 

avi / 

Air Conditioning: 

W = 300 + 
ac 

r l 0. 72 

35 D^ • (No. Seat Rows ) 

Anti-Ice: 

W . = 6.25 • 
ai 

(b/Cos A o/2 )°' 95 


Auxiliary Gear; 

O 1.55 

W =0. 011 (GW x 10 3 ) 

aux 


Furnishings; 

W r = 40 • PAX 1 * 185 + 1625 
fra 


APU (if present): 

W =29.2 

apu 

• PAX 0,7 
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Crew: W = N • W + N • W A (21) 

cr cr crew stew stew 

where 


w 

cr 

= Crew weight kg (lb) 


N 

cr 

= Number in cockpit crew 


W 

= Weight of crew member kg 

(lb) 

crew 

N j 

= Number of stewardesses 


stew 

- 


W , 
stew 

= Weight of stewardesses kg 

(lb). 


The weights of crew members and stewardesses are input values, as is the num- 
ber of flight crew members. The number of stewardesses required is computed 
as a function of the number of passengers in coach and first class sections. In 
first class, there is one stewardess for each 20 passengers; in coach, one for 
each 40 passengers. 


Equipment: W = 100 + 30 • PAX 

eq 

Oil: W = 14.06 ♦ N 

oil eng 

Unusable Fuel: W = 0. 025 * S 

uf w 


( 22 ) 

(23) 

(24) 


In the above three equations 

Passenger equipment weight kg (lb) 

Total number of passenger seats (input) 

Weight of engine oil kg (lb) 

Number of engines (input) 

Weight of scaled engine kg (lb) 

Weight of unusable fuel kg (lb) 

2 2 

Wing area m (ft ) . 

11. Reserve Fuel — FAA regulations define reserve fuel requirements sufficient for 
a holding period at the end of cruise plus diversion to an alternate destination. 
The holding period and distance to alternate vary for domestic and international 
service; these two parameters are therefore defined by input to the model. 


W 

eq 

PAX = 


W 

oil 
N 


eng 


W 

i 

W 


eng 

uf 


w 


2-17 



Reserve fuel for hold at end of cruise is computed by adding the equivalent dis- 
tance to design range. Diversion to an alternate airfield presumably occurs 
following a normal descent; reserve fuel requirements thus include climb to and 
cruise at some intermediate altitude. An altitude of 4573 m (15, 000 ft) is built 
into the model. 


Alternate field reserve fuel requirements are computed in the following manner. 
Since the total fuel supply would be exhausted at the end of this leg, airplane 
gross weight at that point is just OWE plus payload. The climb/descent routine 
described in the following section determines fuel and distance for the descent 
from 4573 m (15, 000 ft), giving gross weighty at the j end of the alternate cruise leg. 
This weight is increased by 5 percent as an estimate of gross weight at the start 
of the diversion. 

Climb fuel and distance to 4573 m (15 , 000 ft) are computed from this weight. 

Climb and descent distances are subtracted from the required alternate field dis- 
tance to yield the necessary cruise distance. 


Range factor is computed at design Mach number at 4573 m (15, 000 ft); the 
method for computing L/D at off-design conditions is described in the following 
section. An inversion of equation (1) is employed to give the ratio of gross 
weights at beginning and end of the alternate cruise leg: 


W /W = exp (D /RF)withD = cruise distance km (nm) . 
o 1 cr cr 

Cruise fuel for this leg is then simply 


W 


fait 


= W (W /W, - 1) 
o 1 


(25) 


where W is the gross weight upon arrival at 4573 m (15,000 feet). Total reserve 
fuel requirement for the alternate field diversion is the sum of Wf^ plus climb 
and descent fuel. This total is added to OWE and payload to give the true gross 
weight at the end of the descent to the primary destination. 

12. Climb and Descent — The method selected for computing climb or descent flight 
path angle, Y, is the simple expression 

Tan Y = T/W - D/L (26) 

which ignores any longitudinal acceleration effects during climb. The same ex- 
pression is valid for computing level flight acceleration during the climb schedule, 
which simplifies programming. 
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The climb/descent schedule employed in the program is shown in Figure 2-7. The 
climb to 3050 m (10, 000 ft) is flown at a constant 130 m/sec (250 knots) equivalent 
airspeed per current federal regulations . The airplane accelerates to maximum 
climb speed, V CJ 0 (input in kt EAS), at 3050 m (10, 000 ft) and climbs at constant 
equivalent airspeed to the altitude at which design Mach number is encountered. 

The remainder of the climb is at constant Mach. 

Gross weight at the start of climb is takeoff weight (GW) less takeoff fuel (W^ Q ), 
which is given by the empirical relation 

\ 

= 0.0152 • N • RATING (27) 

fto eng 

where N en g is the number of engines and RATING is the sea level static thrust of 
the scaled engine. For those instances where water injection is used to reduce 
nitrous oxide omissions during takeoff, the weight of water expended is also sub- 
tracted from GW. This weight is also a function of engine thrust: > 

W TT = 0. 012075 • N • RATING . (28) 

H 2 0 eng 

Climb distance increments, AD, are derived from the relation 

AD = AH/Tan Y ' (29) 

where AH is the altitude increment for a given leg and Tan Y is an average of the 
values computed at the end points of the leg. Elapsed time is given by 

AT = AD/V t (30) 

and the weight of fuel consumed by 

AW f = AT • F • c> (31) 

f n 

where V t is true airspeed, F n is net engine thrust, and c’ is specific fuel con- 
sumption . All three are average values for the leg. 

To compute flight path angle, engine thrust and fuel consumption are interpolated 
from input tables. To obtain D/L, lift coefficient is given by 

C „ = W/q S 
Z w 

where W is the current gross weight and q is dynamic pressure N/m 


(32) 

(lb/ft 2 ) . 
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Figure 2-7. Climb and Descent Schedule 
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At 130 m/sec (250 kt) EAS, q = 10,150 N/m 2 (212 lb/ft 2 ). At 3050 m (10, 000 ft), 
q (lb/ft 2 ) = 3.392 x 10“ 3 V^ 2 . These two values suffice until the constant Mach 
climb segment is reached; thereafter, 

q = 1481 • -6 • M 2 ~ (33) 

with all parameters previously defined. 

Drag coefficient is derived as a function of incremental Mach (AM) below critical 
Mach. Figure 2-8 shows the curves employed. For supercritical airfoils, critic- 
al Mach number is independent of lift coefficient; increments in parasitic and in- 
duced drag coefficients (AC^ and AC^.) are obtained directly from tables repre- 
senting these curves. For conventional 1 airfoils, critical Mach number varies in- 
versely with lift coefficient. For reasonable values of G^, incremental induced 
drag coefficient, AC^. can be considered independent of this effect. To derive 
the increment in parasitic drag coefficient, A Cd 0> the increase in critical Mach 
number (AM C ) for the current is computed. The ACd 0 table is then entered at 
(AM C ) and at (AM C + AM). The difference in these two values interpolated from 
the table is the required AC^ . Total drag coefficient is them 

C J = + AC + C 2 /n M e + AC, (34) 

d a a X d 

o o 1 

Wing aspect ratio, At, and Oswald efficiency factor, e, are input values. 

Returning to the main climb/descent routine, the climb from 457.3 to 3050 meters 
(1500 to 10, 000 feet) is computed as a single increment, first at constant weight 
and again with corrected for the computed change in gross weight. A similar 
technique is employed for the level acceleration leg. The climb at constant V c ^ 
is broken into increments by Mach number since drag coefficient changes rapidly 
during this leg. The Mach number increment used is an input variable, with a 
recommended value of 0. 05. Each increment is computed at constant weight with 
no iteration. The constant Mach climb is likewise broken into incremental alti- 
tude steps because of rapidly changing q. The altitude increment is also under 
external control with a recommended value of 1524 meters (5, 000 feet). 

Descents are computed as climbs, but using idle thrust rather than climb thrust. 
Entering weight for the routine is actually the weight at the end of the descent, 
with fuel being added back in as the descent is computed backwards, starting at 
457.3 m. (1500 ft) up to final cruise altitude. Commercial airplanes normally 
perform cruise step climbs, as fuel is burned off, for reasons of air traffic con- 
trol. The Breguet equation theoretically presumes a continuous climb to main- 
tain constant L/D during cruise. The model follows this assumption and computes 
a final cruise altitude which gives the same cruise L/D at the final cruise weight , 
W 1 , Constant L/D implies constant , and since M is also constant, final cruise 
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altitude can be obtained from the table of altitude vs. Vg/M. At initial cruise 
conditions (subscript o) 

V 2 = 2W /p C„ S . 

e o o SL w 

o 

At end of cruise (subscript 1), only the weight has changed, therefore 

V 2 = 2W /p C S . 

e^ 1 o H w 

By looking up the value of V / M at initial cruise altitude, 

6 


(V /Mi = (V /M) VwVW . 
e 1 e o 1 o 


(35) 


Final cruise altitude is then a function of (V /M), , obtained from the table. 

el 


13. Gross Weight Solution — All of the prerequisite pieces are now available for a 
gross weight solution. The weight ratio W q /W^ is obtained from 


W SW = exp (D cr /RF) 


(36) 


where RF is an average range factor for initial and final cruise conditions and 
D q is cruise distance, given by 


D = R - D . - D J + D 
cr des res 


(37) 


where D c ^ is total climb distance and D(j es is total descent distance. All distance 
computations are in units of nautical miles. The term, D res , is the distance 
equivalent of the reserve fuel loiter time requirement plus a 6 minute allowance 
for air maneuver. The final cruise weight W-^ is known as the weight at the be- 
ginning of the descent. Initial cruise weight is then available from equation (36) 
and gross weight at start of takeoff, GW, is given by 


GW = W + 


W, . + 

fcl 


W + 
fto 


W. 


H 2° 


(38) 


The portion of cruise fuel corresponding to D is computed from 

res 

w fres - W 1 h < D res /RF > ' F ] . 


(39) 


Total fuel capacity is computed by summation of the various increments for climb, 
descent, cruise, and reserves. 
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A single program iteration is performed beginning at step 4 , but omitting step 8 . 


Output — Block fuel and block time computations are performed for output presen- 
tation and for use in the DOC and ROI routines which follow. Block time, Tj^, is 
given in hours by 


T = t + T + T 
bl cl cr des 


+ 0.1 


(40) 


where T c j and T^gg are the computed climb and descent time, the 0. 1 hour con- 
stant is air maneuver time, and cruise time (T cr ) is the quotient of cruise dis- 
tance and average true airspeed. At design range, block fuel, W^fj, is just total 
fuel minus reserves. 


The DOC and ROI routines must determine block time and block fuel as functions 
of distance flown. For this purpose, slope-intercept parameters are computed 
for each of these functions as shown below. 

Block time 

Slope: 

Intercept: 

Block fuel 

Slope: 

Intercept: 

where 

T = Time (hr) 

D = Distance km (nm) 

= Fuel weight kg (lb) 

and subscripts 


T = T /D hrAm (hr/nm) 
bs cr cr 

(41) 

T = T, - R • T, hr 
bi bl bs 

(42) 

F bs = W fcr /D cr kk/km flb/nm) 

(43) 

F bi = W fbl - H ' F bs k S <“» 

(44) 


cr = Cruise 
bl = Block 


A complete sample output is included in Section 5 



2. 2 DOC (DIRECT OPERATING COST) 


The basic function of the DOC module is to translate aircraft characteristics into gen- 
eralized economic dimensions. This is accomplished by a series of equations (cost 
estimating relationships or CERs) which are responsive to design characteristics, 
performance, and operating rules. The first major step is calculation of aircraft fly- 
away cost. 


Flyaway Cost 


Flyaway cost is computed at either of two levels, dependent on availability of design 
definition. The first level computes airframe cost according to the formula 


Cost Airframe = WAF 



98 

0.135 


5500 

PQ 


200,000 

WAF 


0.2 


where 


WAF 

PQ 

98 

0.135 

5500 

0,2 

200, 000 


= Airframe weight kg (lb) 

= Production quantity 

= Nominal $98/pound recurring production cost at unit #1 ($44. 5/kg) 
= Learning curve exponent 

= Nominal $5500/lb development cost ($2500/kg) 

= Cost/size scaling exponents 

= Nominal airframe reference weight in pounds (90909 kg) , 


Engine costs are input directly or computed as a function of thrust. Avionics costs are 
input directly. Total flyaway cost is the sum of airframe, engines, and avionics costs. 


When functional weights and cost coefficients are available, airframe costs are com- 
puted in 24 functional elements by labor and material as follows: 


Cost of Airframe Element = Labor Cost + Material Cost 


Labor Cost = 



MLL 


CLSAL 

• WSUB • C MLH * (CRCS) • CLRHW 


where 

PQ = Production quantity actual 

PR = Production quantity reference 

MLL = Labor learning exponent 
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WSUB = Subsystem weight kg (lb) 

CMLH = Cost per labor Hour — manufacturing $/hr 
CRCS = Reference weight kg (lb) 

CLSAL = Labor size scaling exponent 

CLRHW = Hours per pound at reference weight I hr/kg (hr/lb) 


Cost of Mat’l 


where 



MML 


• WSUB • MCF • (CRCsf MSAL • CMRCW 


MML = Material quantity learning factor 

MCF = Material unit cost at reference weight $/kg ($/lb) 

CMSAL = Material size scaling exponent 

CMRCW = Material unit cost at reference weight $/kg ($/lb) 

MLL — Labor Learning Exponent 

In the labor cost equation, the exponent MLL is used to adjust from one production 
quantity to another along a cost reduction or learning curve. This approximation is 
intended for use between 100 units and 2000 units, and is not as accurate as rigorous 
learning curves applied back to the first production unit. Figure 2.2-1 shows expon- 
ents approximating various learning curves from unit 100 to unit 500 — the principal 
region of interest. For the general case, a labor curve of 80% is recommended. 

This corresponds to a value at MLL = 0.680. 

CLSAL — Labor Sizing Scaling Exponent 

For either a given component (e. g. , a wing) or for a total airframe, there is signifi- 
cant cost/size scaling. That is, as design size is increased, cost per unit of weight 
is decreased. Figure 2.2-2 shows transport cost vs. weight. Figure 2.2-3 shows 
horizontal stabilizer hours as a function of weight. In these two examples, the size 
scaling exponents are both 0.8. 

For labor only, the following size scaling exponents have been developed on prior 
programs: 


1. Total Airframe 0. 82 

2. Primary Structure — Wing, Fuselage, Empennage 0,80 

3. Electrical, Pneumatic, Hydraulic 0.85 

4. Final Assembly and Checkout 0,70 
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MML — Material Quantity Learning Factor 

Provision is made in the model for use of a material quantity cost reduction curve 
where the analyst has cost/quantity data for material. This "lear ning 11 is much less 
than labor cost/quantity reduction and in the absence of detailed data, a 98% curve 
(MML = 0.96) is suggested as being typical. 

CM SAL — Material Size Scaling Exponent 

This factor is available for use where detailed cost data is available on material as a 
function of airframe or component size. With typical, light alloy construction, this 
factor is very close to 1.0; however, with composite materials, special skin gages, 
or specially tailored alloys, scaling may be required. For some applications requir- 
ing highly specialized machinery or if a critical fraction of some resource is required, 
it might take on a value greater than one. 

4 

sustaining tooling and engineering are computed as percentages of initial tooling and 
engineering. Assembly, integration, and profit are computed as percentages. Initial 
tooling and engineering are pro-rated and added to establish total airframe cost. The 
sample run and listing contains the equations and a set of input values representative 
of a conventional large transport aircraft. 

Direct Operating Costs 

Direct operating costs are computed according to the 1967 ATA formula, with cost 
factors updated to 1971 rates. Cost factors requiring updating for the effects of infla- 
tion are as follows; 

CREW = Crew cost constant $/hr — 1971 Value $180 
CDLH = Maintenance direct labor $/hr - 1971 Value $5 

CFT = Fuel cost $/metric ton — 1971 Value $32. 60/metric ton ($15/1000 lb) 

Costs are computed for the various distance specified in DNM(I) — up to 20 separate 
distances. Times are computed from parasitic time (taxi, takeoff, air maneuver) 
identified as TBI (block time intercept from synthesis) and TBS (block time slope 
from synthesis). Block fuel is obtained from the synthesis and adjusted for fuel used 
in air delay (AD) and ground delay (GD), if any. Speeds, times, and block fuel are 
calculated for each DNM(I). 

Aircraft utilization is calculated in an approximation of the ATA formula as follows; 

U = ANHR • TBTR/ (TURNT + TBTR) 
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where 


U = Annual utilization block hours 
ANHR = Idealized utilization (4800 hr/yr) 

TBTR = Block time for a given flight hr 

TURNT= Ground turnaround time at gate (typically 0.50 hours) 

Figure 2-9 shows resulting annual block utilization. 

The ATA equations are then used to calculate Direct Operating Costs expressed in 
dollars per nautical mile for the following elements (see listing for equations). 

CCT = Crew cost 

CFOT = Cost of fuel and oil 

CINST = Cost of (hull) insurance 

CAFLT= Cost of airframe and misc. labor (maintenance) 

CAFMT=Cost of airframe and misc. material (maintenance) 

CCELT = Cost of engine labor (maintenance) 

CCEMT- Cost of engine material (maintenance) 

CDEPT = Cost of depreciation, including depreciation of spare parts 
CAMBT= Cost of applied maintenance burden 

These elements are then printed out for each DNM(I) as shown on the sample format. 

2.3 ROI MODULE (RETURN-ON-INVESTMENT) 

Direct operating cost is an abstract measure of an airplane's economic efficiency fly- 
ing an arbitrary distance. ROI attempts to relate economic efficiency to actual city 
pairs, actual traffic levels, and reasonable fare yields. In addition to the direct costs, 
the indirect costs of passenger services, aircraft handling, and administration are in- 
cluded in the calculations. 

The basic returo-on-investment percentage is defined as that interest rate on a fixed 
dollar investment (e. g. , a bank savings account) which will yield equal financial bene- 
fit over the investment period to the investment being analyzed. 

For an airplane, the investment period is the depreciable lifetime. The formula is 

ROIF = R • (1 + R) Y ^*/'((1 + R)™ 5 -!) 
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Annual 

Block 

Utilization 

(Hours) 



ROIF = ROI percentage 

R = Annual percentage of profit on investment 
YRS = Years of investment term . 

An iterative solution is done by computer yielding an answer within a prespecified 
error specified by EPSA (typically 0. 001). The iteration proceeds by a stepping factor 
in the iteration specified by DELR, typically 0. 05. 

To obtain R, the annual percentage of profit, it is necessary to calculate the following: 

1. Flights per day by city pair to handle input traffic 

2. DOC by city pair 

3. IOC by city pair 

4. Income by city pair 

5. Net income after taxes 

6. Aircraft required and total investment 

Flights per day per city pair are computed by establishing a minimum frequency, a 
maximum permissible average load factor, and dividing traffic (input) by capacity. 
Number of flights are rounded to next higher integer, giving resulting load factors be- 
low maximum specified load factor. Referring to the sample printout, it is seen that 
low traffic city pairs are constrained by the minimum frequency and have low load 
factors. City pairs with greater activity, say 10 or 20 flights per day approach the 
maximum specified load factor — in this case 55%. 

DOCs are computed for the specified city pair distances by interpolation from DOCs 
at arbitrary distances calculated in the previous section. IOCs are calculated as 


follows, based on the Lockheed IOC methodology. 

Ramp Service and Landing Fee $3. 02 per metric ton 

($1.37 per 1000 lb) gross weight per landing Z (1) 

Passenger Service Cost — $0.0047 per RPM Z (2) 

Reservations, Sales, and Passenger Handling — $9/passenger 

hour Z (3) 

Stewardess Expense — $15. 03/Stewardess Block Hour Z (4) 

(No. of Stewardesses from Synthesis) 

Meal Cost = $0.90 + $ 0.45/Passenger Block Hour Z (5) 
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z (6) 


Administrative Cost Factor = 0. 091 


IOC Costs 

1 , 


All of the above are in 1971 dollars and are subject to inflation, except the adminis- 
trative override. 


Traffic data for each of up to 200 city pair is entered as: 


City Pair 
Distance 
Pass. /Day 


See input format 


Income by city pair is calculated from a fare formula as follows: 


For each passenger $7 (FARE I) 

Plus $0. 0346/km ($0. 0644/nm) (FARES) 

This was derived empirically for 1971 by talcing actual trunk airline income and RPM, 
and as such, accounts for cargo, mix of first class and coach passengers, and dilutions 
in the form of special and promotional fares. Fares must be adjusted for inflation, 
and the user is cautioned to look at actual yields, not published fare schedules, to 
develop realistic income levels. Typical yields are about 85% of published regular 
coach fares. 

Net income after taxes is computed by taking total income less DOC, less IOC times 52%. 
This is then defined as net profit. The total aircraft required are summed from frac- 
tional requirements by city pair. 

This is done on the assumption that schedules can be arranged to utilize airplanes 
within the constraints of the empirical utilization formula, which was presumably 
developed from actual data. For each aircraft required, investment is equal to fly- 
away cost plus engine and airframe spares. Total investment is obtained by multiply- 
ing investment base per airplane (XLB) by fleet size. Total investment is compared to 
system net profit to obtain system Retum-on-Ihvestment. 

Table 2.2 presents a summary of terms used in the analytic derivation described in 
Section 2 of this report. Both the scientific international and the English system of 
units are presented; however, the actual program utilizes the English system of units 
' exclusively. ' ~ 
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Table 2. 2 Summary of Terms 


A 


A 

cs 




A 

nac 

ANHR 



AR 


b 


c 

F 

c 


CAFLT 

CAFMT 

CAMBT 

CCELT 

CCEMT 

CCT 

CDEPT 

CDLH 


Total airplane wetted area m (ft ) 

2 2 

Airfoil cross sectional area m (ft ) 

2 2 

Fuselage gross wetted area m (ft ) 

2 2 

Wing inboard portion wetted area m (ft ) 

2 2 

Nacelle wetted area m (ft ) 

Idealized utilization 4800 hr/yr 

2 2 

Wing outboard portion wetted area m (ft ) 

Aspect ratio (dimensionless) 

Wing span m (ft) 

Wing chord m (ft) 

Specific fuel consumption kg/N-hr (lbm/lb-hr) 
Cost of airframe maintenance labor $ 

Cost of airframe maintenance material $ 

Cost of applied maintenance burden $ 

Cost of engine maintenanc e labor $ 

\ 

Cost of engine maintenance material $ 

Crew cost $ 

Cost of depreciation including spares $ 
Maintenance direct cost per hour $/hr 
Equivalent skin friction coefficient (dimensionless) 
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CFOT 

CFT 


cgl 

CINST 

°i 

CLRHW 

CLSAL 

V 

c di 

C do 

CMLH 

CMRCW 

CMSAL 

CRCS 

CREW 

C . 
root 

c 

sw 

D 

D 

0 
D 

1 

D 


c SL 


cr 


des 


Cost of fuel and oil $ 

Fuel cost $/ metric ton ($/1000 lb. ) 

Ratio of glove root chord to theroretical wing chord (dimensionless) 
Cost of (hull) insurance $ 

Lift coefficient (dimensionless) 

Hours per pound at reference weight 
Labor size scaling exponent (dimensionless) 

Total drag coefficient (dimensionless) 

Induced drag coefficient (dimensionless) 

Parasetic drag coefficient (dimensionless) 

Manufacturing cost per labor hour $/hr 

/ 

Material unit cost of reference weight $/kg ($/lb.) 

Material size scaling exponent (dimensionless) 

Reference weight kg (lb) 

Crew cost per hour $/hr 
Wing root chord m(ft) 

Ratio wing wetted area to planform area (dimensionless) 

Scaled engine nacelle diameter m (ft) 

Nacelle diameter of 1 . 0 engine m (ft) 

Total- climb distance km-(nm) - - ----- 

Total cruise distance km (nm) 

Total descent distance km (nm) 
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D. 

fus 

DOC 

D 

nac 

D 

res 

e 

EAS 

ES 


Fuselage external diameter m (ft) 

Direct Operating Cost $ 

Nacelle diameter m (ft) 

Total reserve distance km (nm) 

Oswald efficiency factor (dimensionless) 
Equivalent airspeed km/hr (kf) 

Engine scale factor (dimensionless) 


br 

F 

n 

GW 


K 


acs 
K 

ar 

KEAS 


K 


matl 


K 


tech 
L/D 

(L/D) 


nac 


(L/D) 

nac-fxn 


Net engine thrust N (lb. ) 

Gross weight kg (lb) 

Active control system factor (dimensionless) 

Fuselage area ruling factor (dimensionless) 

Knots equivalent airspeed (kt) 

Construction material factor (dimensionless) 

Aerodynamic and structural technology factor (dimensionless) 
Optimum lift to drag ratio (dimensionless) 

Nacelle length to diameter ratio (dimensionless) 

Fin mounted nacelle length to diameter ratio (dimensionless) 



m. a. c. 


Length of fuselage straight section m (ft) 

Length of fuselage tapered nose and tail section m (ft) 

Mean aerodynamic chord m (ft) 

Freestream Mach number (dimensionless) 
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M 



M 

c 

MCF 

MLL 

MML 

N 

cr 

N 

eng 

N * 
stew 

OWE 

PAX 

PQ 

PR 

g 

R 

R 

RATING 

RF 

ROI 

ROIF 

RPM 



S 

v 


Critical Mach number (dimensionless) 

Material unit cost at reference weight $/kg ($/lb) 

Labor learning exponent (dimensionless) 

Material quantity learning factor (dimensionless) 

Number of cockpit crew (dimensionless) 

Number of engines (dimensionless) 

Number of stewardesses (dimensionless) 

Basic operating (empty) weight kg (lb.) 

Total number of passenger seats (dimensionless) 

Actual production quantity (dimensionless) 

Reference production quantity (dimensionless) 

2 2 

Dynamic pressure N/m (lb. /ft ) 

Annual percentage of profit on investment (dimensionless) 
Range km (nm) 

Sea level static thrust of scaled engine N (lb. ) 

Range factor km (nm) 

Return on investment (dimensionless) 

Return on investment percentage (dimensionless) 

Revenue passenger miles km (nm) 

2 2 

Horizontal tail exposed area m (ft ) 

2 2 

Wing inboard portion planform area m (ft ) 

2 2 

Wing outboard portion planform area m (ft ) 

2 2 

Vertical tail exposed area m (ft ) 


S 

w 


2 2 

Wing area m (ft ) 
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TBTR Block time for a given flight (hr.) 

t/ c Airfoil thickness to chord ratio (dimensionless) 

t Wing thickness at the root m (ft) 

TURNT Ground turnaround time at gate (hr) 

T/W Cruise thrust to weight ratio (dimensionless) 

u Glove root chord m (ft) 

U Annual utilization (hr) 

v Basic wing root chord m (ft) 

Equivalent airspeed km/hr (kt) 

V, Horizontal tail volume coefficient (dimensionless) 

h 

V Time airspeed km/hr (kt) 

V 

V Vertical tail volume coefficient (dimensionless) 


W 

W 

o 


W 

ac 

WAF 

W . 
ai 

W 

apu 

W 

aux 

W . 
avx 

WAVI 


Current gross weight kg (lb.) 

Initial cruise gross weight kg (lb. ) 

Final cruise gross weight kg (lb. ) 

Air conditioning system weight kg (lb . ) 
Airframe weight kg (lb. ) 

Anti-icing system weight kg (lb. ) 

Auxiliary power unit system weight kg (lb. ) 
Auxiliary gear weight kg (lb.) 

Avionics system weight kg (lb. ) 

Weight of avionics kg (lb.) 
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w. 


w 


bfl 


W 


cr 


W 


crew 


W 


dist 


W_ 

dmp 

WE 


W 


ec 


W 


elec 


W 


eng 


W 


engl 


W 


eg 


W 


fait 


W 


fc 


W 


fcl 


W 


fcr 


W. 


fres 


W 


fra 


W 


fto 


W 


flot 


W, 


W 


h20 


W 


hyd 


Fuselage structural weight kg (lb.) 

Block fuel weight kg (lb. ) 

Crew weight kg (lb. ) 

Weight of a crew member kg (lb. ) 

Fuel system distribution system weight kg (lb. ) 

Fuel system dumping system weight kg (lb. ) 

Weight empty kg (lb. ) 

Engine controls weight kg (lb. ) 

Electrical system weight kg (lb. ) 

Individual engine weight kg (lb.) 

Individual scaled engine weight kg (lb. ) 

( 

Passenger equipment weight kg (lb. ) 

Cruise fuel weight for alternate cruise leg kg (lb.) 
Fuel system controls weight kg (lb. ) 

Climb fuel weight kg (lb. ) 

Cruise fuel weight kg (lb.) 

Fuel weight corresponding to reserve distance kg (lb. ) 
Furhishings weight kg (lb. ) 

Takeoff fuel weight kg (lb.) 

Total fuel weight kg (lb. ) 

Horizontal tail structural weight kg (lb. ) 

Water injection system fluid weight kg (lb.) 

Hydraulic system weight kg (lb. ) 
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w 


inst 


W 


Jpg 


W 

oil 

W \ 
nac 

W 


W 


pl 

W _ 
rf 

W/S 


w 


sc 


W 


si 


W 


sp 


W 


st 
W 

stew 

WSUB 


W 


uf 


W 


W 


W 


w 


w . 

W1 


XIB 


gl 

YRS 


Instrument weight kg (lb. ) 

Landing gear weight kg (lb.) 

Engine oil weight kg (lb. ) 

Nacelle weight kg (lb. ) 

Fuel system pump weight kg (lb. ) 

Payload weight kg (lb.) 

Fuel system refueling system weight kg (lb. ) 

2 2 

Design takeoff wing loading N/ m (lb/ft ) 

Surface control system weight kg (lb. ) 

Fuel system sealant weight kg (lb. ) 

Engine/ nacelle soundproofing weight kg (lb. ) 
Engine starting system weight kg (lb. ) 
Stewardesses weight kg (lb. ) 

Subsystem weight kg (lb.) 

Unusable fuel weight kg (lb.) 

Fuel system vent system weight kg (lb. ) 

Vertical tail structural weight kg (lb. ) 

Wing structural weight kg (lb. ) 

Engine water injection system weight kg (lb. ) 
Investment base per airplane $ 

Location of wing root/fuselage intersection nx (ft.) 
Years of investment term (yr. ) 


2-38 



y 

6 

AD 


AH 

AT 

A c/2 

^LE 


Climb or descent flight path angle (degrees) 

Atmospheric pressure ratio at initial cruise altitude (dimensionless) 
Climb distance increment km (nm) 

Climb altitutde increment m (ft.) 

Climb time increment (hr) 

Wing 50% chord sweep (degrees) 

Wing leading edge sweep (degrees) 

3 3 

Atmospheric density at initial cruise altitude kg/m (lb/ft ) 
Atmospheric pressure ratio at a given altitude (dimensionless) 
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SECTION 3 


COMPUTER PROGRAM 


3.1 PROGRAM ORGANIZATION 

The TEKVAL computer program is configured for maximum flexibility of operation. 
Each of the main program sections (SIZE, DOC, and ROI) may be operated indepen- 
dently with external input, or they may be operated sequentially (SIZE-DOC-ROI or 
DOC -ROI) with the output of one section providing partial input for subsequent sections. 

Program operation is controlled by a 3-digit binary control word labeled KEY. 

Figure 3-1 illustrates the TEKVAL program organization and the function of the 
control word. The three digits control the execution of each of the main program 
modules, the first digit for SIZE, the second for DOC, and the third for ROI. A one 
calls for module execution while a zero signals omission of that module. To operate 
all three modules sequentially, a value of KEY = 111 is used. A value such as KEY = 
010 calls for execution of the DOC module only. 

The value of KEY is read at the beginning of the main program. A branch is then 
made to the appropriate module within the main program. Following execution of 
each module, the value of KEY is reduced by an appropriate power of 10, effectively 
eliminating the leading digit. If the succeeding digit so indicates, the next module is 
executed. If the succeeding digit is a zero, the program reverts to the beginning for 
a new value of KEY. 

Each of the main program modules has its own input and output subroutines. The 
NAMELIST format is utilized for all input, including the control word KEY. A 
complete input deck is obviously required for execution of the first module of a given 
run. For sequential operation of subsequent modules, only partial input decks are 
required but there must be a separate NAMELIST deck for each module to be executed 
on every run. Complete input requirements and data communication between modules 
are described below. 

3. 2 DESCRIPTION OF SUBROUTINES 

Figures 3-2 through 3-4 show the subroutines called by each of the main program 
modules and how . they are interlinked. All input and output subroutines are straight- 
forward and need no further clarification. The principle function of each of the 
remaining subroutines is described herein. 
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Figure 3-1. TEKVAL Program Organization 
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Figure 3-2. SIZE Subroutine Structure 
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Figure 3-3. DOC Subroutine Structure 













The following subroutines are employed by the SIZE module: 

TLUI - A one dimensional linear interpolation table look-up routine. The call 
sequence lists the input value X and the locations and dimension of the X and Y 
vectors. The value of Y corresponding to X is returned via a call argument. If 
the vector dimension is exceeded during the table search, an error message "TLUI 
OUT OF RANGE" is printed and the program stops. 

TLU2 - A two dimensional linear interpolation table look-up routine. This routine 
gives Y = f (X, Z). The call sequence lists values of X and Z plus addresses for the 
X and Z vectors and the Y matrix, and the dimensions of the X and Z vectors. The 
Z indices are determined first. For each of the two Z indices, the TLUI subroutine 
is utilized to determine Y = f(X). Linear interpolation on Z then gives the final 
value of Y. Should the Z vector dimension be exceeded, the error message "TLU2 
OUT OF RANGE" is printed and the program stops. 

CLDES - This subroutine computes the time, distance, and fuel required to climb or 
descend according to a given speed-altitude profile. The mathematical operations 
were described in Section 2.1. This routine performs the computations and accumu- 
lations of distance, time, and fuel increments as well as current gross weight, Mach 
number, dynamic pressure, true airspeed,and average values for climb gradient, 
specific fuel consumption, and engine thrust. It directly utilizes TLUI for altitude 
related functions and frequently calls TANGAM to compute climb gradients. 

TANGAM - This subroutine computes flight path gradient as the difference between 
T/W and D/l. Subroutine DOVRL is called to produce D/L. Subroutine TLU2 is 
called to obtain engine thrust and specific fuel consumption, which are then returned 
to CLDES via TANGAM. 

DOVRL - Subroutine DOVRL computes the ratio C^/C^ as a function of Mach number 
for a given value of C^. The primary function of this subroutine is to back out com- 
pressibility effects from the cruise Mach number drag coefficient for the slower 
speeds encountered during climb and descent. The mathematical approach was 
described in Section 2.1. 

The DOC and ROI modules each employ only one subroutine in addition to the input 
and output functions. Subroutine GENEQ makes use of a generalized cost equation, 
explained in Section 2.2, to compute labor and material costs, and their total, for 
a particular aircraft component or subsystem. The call arguments transfer the 
component weight, a reference weight, labor and material scaling exponents, and labor 
hours/pound and material cost/pound at the reference weight. All of these factors 
vary by component. A COMMON statement provides constant values for desired and 
reference production quantities, labor and material learning exponents, manufacturing 
labor rate, and material cost escalation factor. 
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In the ROI module, subroutine RET employs an iterative technique to compute the 
equivalent rate of return-on-investment equaling the profit generated by commercial 
airplane operations. The call sequence provides all required input: profit, number 
of airplanes,, investment/airplane, and the time period of the operation. The mathe- 
matical technique is described in Section 2.3; because of the nature of the equations, 
an iterative technique is required. A trial value of ROI is used to compute what is 
in essence an equivalent profit for a fixed investment over a given number of years. 
The trial ROI is incremented in steps of DELR (arbitrary stepping increment, dimen- 
sionless) until the computed profit matches the target value within an accuracy 
specified by the parameter EPSA ( error fraction, dimensionless). Both DELR and 
EPSA are input; suggested values are DELR = .05, EPSA = .001. When profit is 
negative, ROI is computed based on the absolute value of profit, with the negative 
sign attached to the result. 

3.3 DICTIONARY OF VARIABLES 

A dictionary of FORTRAN variables found in the program modules associated with the 
sizing process is presented in Table 3-1. Variables associated with the Direct Oper- 
ating Cost and Return-On-Investment modules are found in Tables 3-2 and 3-3, re- 
spectively. It should be noted that the actual program input and output utilize the 
English system of units. 


i 

I 


I 


I 

I 
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TABLE 3-1. DICTIONARY OF VARIABLES 
NAME DESCRIPTION 

SIZE MODULE 


ALT INITIAL CRUISE ALTITUDE (FT/1000) 

ALTF FINAL CRUISE ALTITUDE (FT/1000) 

CDF EQUIVALENT SKIN FRICTION COEFFICIENT 

CDOMCR CRUISE PARASITE DRAG COEFFICIENT 

CGLOVE RATIO WING GLOVE CHORD/THEO CHORD AT BODY 
CHVOL VOLUME COEFFICIENT FOR HORIZONTAL TAIL 
CL LIFT COEFFICIENT 

CLCR CRUISE LIFT COEFFICIENT 

CLOCH HORIZONTAL TAIL MAC LOCATION RATIO TO BODY LENGTH 
CLOCV VERTICAL TAIL MAC LOCATION RATIO TO BODY LENGTH 
CLOCW WING MAC LOCATION RATIO TO BODY LENGTH 

CMDD INCREMENT IN CRITICAL MACH NUMBER DUE TO CL / CL 

CPAXC FRACTION OF TOTAL SEATS IN COACH 

CRMACH DESIGN CRUISE MACH NUMBER 

CSWET RATIO WING WETTED AREA/PLANFORM AREA - RECTANGULAR SECTION 

CSWPYL RATIO PYLON/NACELLE WETTED AREAS 

CTWR INITIAL ESTIMATE OF STATIC T/W REQUIRED 

CVVOL VOLUME COEFFICIENT FOR VERTICAL TAIL 

CWEF INDUCED DRAG EFFICIENCY FACTOR 

CWFRG LATERAL EXTENT OF WING GLOVE/WING SEMI-SPAN 

CWT RATIO WO/WGTO 

CWTIC RATIO WO/WGTO .INITIAL ESTIMATE) 

DAISLE AISLE WIDTH IN COACH (IN) 

DBD BODY DIAMETER (FT) 

DBL TOTAL BODY LENGTH (FT) 

DBLC LENGTH OF COACH SECTION (FT) 

DBLST LENGTH OF STRAIGHT BODY SECTION (FT) 

D8LTPR LENGTH OF TAPERED BODY SECTION (FT) 

DBL1 LENGTH OF FIRST CLASS SECTION (FT) 

DCDCV TABLE - DRAG RISE - DELTA CDO - CONVENTIONAL (COUNTS) 

DCDICV TABLE - DRAG RISE - DELTA CDI - CONVENTIONAL (COUNTS) , 

DCDISC TABLE - DRAG RISE - DELTA CDI - SUPERCRITICAL (COUNTS) 

DCDSC TABLE - DRAG RISE - DELTA CDO - SUPERCRITICAL (COUNTS) 

DCL DISTANCE - CLIMB TO INITIAL CRUISE ALTITUDE (NM) 

DCR CRUISE DISTANCE ( NM ) 

DCRES CLIMB DISTANCE TO 15000 FT FOR FUEL RESERVES (NM) 

DDES DISTANCE - DESCEND FROM FINAL CRUISE ALT ITUOE (NM) 

DDRES DESCENT DISTANCE FROM 15000 FT FOR FUEL RESERVES (NM) 

DRES2 CRUISE DISTANCE TO ALTERNATE DESTINATION (NM) 

DELTAI PRESSURE RATIO AT INITIAL CRUISE ALTITUDE 

DHL LENGTH OF HORIZONTAL TAIL. ARM (FTL 

DLCR CRUISE D/L RATIO 

DLNAC LENGTH OF NACELLE (FT) 

DMACH TABLE - DRAG RISE - MACH NUMBER 

DNAC NACELLE DIAMETER(FT) 

DNAC1 DIAMETER OF FULL SCALE NACELLE (FT) . 

DRES DISTANCE TO ALTERNATE DESTINATION (NM). 

DRESl DISTANCE EQUIVALENT OF TRES (NM) 

DSTPC SEAT PITCH IN COACH (IN) 

DSTP1 SEAT PITCH IN FIRST CLASS (IN) 

DSTWC SEAT WIDTH IN COACH (IN) 

DSTW1 SEAT WIDTH IN COACH (IN) 


DWB WING SPAN (FT) 
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NAME 

owes 

DWCGL 

DWCOB 

DWCR 

DWCT 

DWMAC 

DWTGL 

DWTR 

DVL 

ES 

FBI 

FBS 

FGAM 

FGAMDD 

fldens 

FNCL 

FNCR 

fnidl 

fnmcr 

FNREQ 

FNRTG 

FNZL 

FTAPER 

GNFR 

GNFRF 

GWAR 

GWPL 

GWS 

GWSWPO 

GWSWP2 

GWSWP4 

GWTC 

GWTPR 

H 

HDEL 

HH 

HSIGSR 

KACS 

KAF 

KAPU 

KARULE 

KEF IN 

KEY 

kmatl 

KPSP 

KRPT 

KTEST 

KUPDN 

KWI 

NAISLE 

NCREW 

NENG 

NENGPR 

NPAXC 

npaxi 


TABLE 3-1. DICTIONARY OF VARIABLES (CONT) 
DESCRIPTION 

WING STRUCTURAL SPAN (FT) 

LENGTH OF WING GLOVE CHORD (FT) 

LENGTH OF WING CHORD AT BODY JUNCTION (FT) 

LENGTH OF WING ROOT CHORD (FT) 

LENGTH OF WING TIP CHORD (FT) 

LENGTH OF WING MEAN AERODYNAMIC CHORD (FT) 

THICKNESS OF WING GLOVE CHORD (FT) 

WING ROOT CHORD THICKNESS AT C/L (FT) 

LENGTH OF VERTICAL TAIL ARM (FT) 

ENGINE SCALE 

BLOCK FUEL EQUATION - INTERCEPT 

BLOCK FUEL EQUATION - SLOPE 

TABLE - CRITICAL MACH INCREMENT - FACTOR 

FACTOR TO DETERMINE CMDD = F(WING SWEEP) 

FUEL DENSITY (LB/GAL) 

TABLE - NET THRUST - CLIMB 
TABLE - NET THRUST - CRUISE 
TABLE - NET THRUST - IDLE - 

FULL SCALE ENGINE NET THRUST - INITIAL CRUISE (LBS) 
REQUIRED NET THRUST/ENGINE - INITIAL CRUISE (LBS) 

SLST OF FULL SCALE ENGINE 

LIMIT LOAD FACTOR 

TAPERED BODY LENgTH(DIAMETERS ) 

PODDED NACELLE FINENESS RATIO 

nacelle fineness ratio - fin mounted 

WING ASPECT RATIO 

RATIO WGTO/WPL 

WING LOADING (LB/SQ FT) 

WING SWEEP ANGLE AT LE (DEG) 

WING SWEEP ANGLE AT C/2 (DEG) 

WING SWEEP ANGLE AT C/4 (DEG) 

WING THICKNESS/CHORD RATIO 
WING TAPER RATIO 

TABLE - STANDARD ATMOSPHERE - ALTITUDE ( FT/1000 > 

TABLE - STANDARD ATMOSPHERE - DELTA 
TABLE - ENGINE DATA - ALTITUDE (FT/ 1000) 

TABLE - STANDARD ATMOSPHERE - SQ RT SIGMA 
FLAG FOR ACTIVE CONTROL SYSTEM (1*WITH, Q=W/0) 

FLAG FOR SUPERCRITICAL AIRFOIL (1=WITH* 0*W/0> 

FLAG FOR AUXILLIARY POWER UNIT (1=WITH* 0*W/0) 

FLAG FOR BODY AREA RULING ( 1 = WI TH * 0= W/O ) 

FLAG FOR FIN MOUNTED ENGINE (1«YES* 0«NO) 

PROGRAM CONTROL WORD 

FLAG FOR CONSTRUCTION MATERIAL (1=C0MP* 0*ALUM) 

FLAG FOR NACELLE SOUND PROOFING (1»WITH* 0=W/0) 

FLAG FOR SIZE ITERATION LOOP 

TEST FOR ODD NUMBER OF ENGINES (+=EVEN, -=ODD) 

FLAG FOR CLDES SUBROUTINE (1=CLIMB* 0=DESC£ND) 

FLAG FOR ENGINE WATER INJECTION PROVISIONS (1=WJTH, 0=W/0) 
NO. AISLES 

NO. OF FLIGHT CREW MEMBERS 
NUMBER OF engines 

NO. OF PODDED ENGINES (WING OR BODY MOUNTED) 

NO. SEATS IN COACH 
NO. SEATS IN FIRST CLASS 
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TABLE 3-1. DICTIONARY OF VARIABLES (CONT) 


NAME 

NROWC 

NROW1 

NSROWC 

NSROW1 

NSTEW 

OWE 

PAX 

PIAE 

Q 

QFUEL 

RACS 

RANGE 

RANGEF 

rating 

RF 

RFAVG 

RMACH 

RMATL 

SBW 

SBWCS 

SFC 

SFCIDL 

SFCCL 

SFCCR 

SHEXP 

SNW 

SRSIG 

STW 

SVEXP 

SW 

SWET 

SWP 

SWW 

SWWIB 

SWWOB 

TBI 

TBL 

TBS 

TCL 

TCR 

TDES 

TGAMC2 

TRES 

TSTO 

VCL 

VDIVM 

VMO 

VM1 

VOVRM 

VT 

VTAVG 
VTF 
VT I 
WAF 
WCRES 


DESCRIPTION 

NO. SEAT ROWS IN COACH 

NO. SEAT ROWS IN FIRST CLASS 

NO. SEATS ABREAST IN COACH 

NO. SEATS ABREAST IN FIRST CLASS 

NO. OF STEWARDESSES 

OPERATING WEIGHT EMPTY (LBS) 

TOTAL NO. SEATS 

PI X ASPECT RATIO X INDUCED DRAG EFF FACTOR 
DYNAMIC PRESSURE (LB/SQ FT) 

FUEL QUANTITY - TOTAL (GAL) 

ACS CORRECTION FACTOR FOE WING WEIGHT 
DESIGN RANGE ( NM ) 

BREGUET RANGE FACTOR 

SEA LEVEL STATIC THRUST OF SCALED ENGINE (LBS) 

RANGE FACTOR AT 15000 FT CRUISE OR AT FINAL CRUISE 
RANGE FACTOR - AVERAGE INITIAL AND FINAL CRUISE 
SPEED CORRECTION FACTOR FOR WING WEIGHT 
MATERIAL CORRECTION FACTOR FOR WING WEIGHT 
WETTED AREA - BODY ( SQ FT) 

WING GLOVE CROSS SECTION AREA ( SQ FT) 

THRUST SPECIFIC FUEL CONSUMPTION (LB/LB/HR) 

TABLE - SPECIFIC FUEL CONSUMPTION - IDLE 
TABLE - SPECIFIC FUEL CONSUMPTION - CLIMB 
TABLE - SPECIFIC FUEL CONSUMPTION - CRUISE 
EXPOSED AREA OF HORIZONTAL TAIL ( SQ FT) 

WETTED AREA - NACELLES (SQ FT) 

SQUARE ROOT OF DENSITY RATIO * F ( ALT ITUDE ) 

WETTED AREA - HORIZONTAL + VERTICAL TAIL (SQ FT) 

EXPOSED AREA OF VERTICAL TAIL ( SQ FT) 

WING AREA (SQ FT) 

WETTED AREA - TOTAL AIRPLANE ( SQ FT) 

TABLE - CRITICAL MACH INCREMENT - WING SWEEP 
WETTED AREA - WING TOTAL ( SQ FT) 

WETTED AREA - WING 2NBOARD SECTION (SQ FT) 

WETTED AREA - WING OUTBOARD SECTION (SQ FT) 

BLOCK TIME EQUATION - INTERCEPT 
BLOCK TIME (HR) 

BLOCK TIME EQUATION - SLOPE 

TIME - CLIMB TO INITIAL CRUISE ALTITUDE (HR) 

TIME - CRUISE (HR) 

TIME - DESCEND FROM FINAL CRUISE ALTITUDE (HR) 

TAN (WING SWEEP AT C/2) 

RESERVE FUEL LOITER TIME (HR) 

ENGINE SLST REQUIRED 

MAX CLIMB EQUIVALENT AIRSPEED (KT) 

TABLE - STANDARD ATMOSPHERE - VE/M 
RATIO VE/M AT START OF CRUISE 
RATIO VE/M AT END OF CRUISE 

RATIO EQUIVALENT AIRSPEED/MACH NUMBER = F( ALTITUDE) 
WORKING VALUE - TRUE AIRSPEED (KT) 

TRUE AIRSPEED - AVG INITIAL AND FINAL CRUISE (KT) 

TRUE AIRSPEED - FINAL CRUISE (KT) 

TRUE AIRSPEED - INITIAL CRUISE (KT) 

AIRFRAME WEIGHT FOR DOC PROGRAM (LBS) 

RESERVE FUEL WEIGHT - ALTERNATE DESTINATION - CLIMB (LBS) 
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TABLE 3-1. DICTIONARY OF VARIABLES (CONT) 


NAME DESCRIPTION 

WCREW WEIGHT OF FLIGHT CREW MEMBER (LBS) 

WDRES RESERVE FUEL WEIGHT - ALTERNATE DESTINATION - DESCENT (LBS) 
WE WEIGHT EMPTY (LBS) 

WENG1 FULL SCALE ENGINE WEIGHT (LBS) 

WFBL BLOCK FUEL WEIGHT (LBS) 

WF CL FUEL WEIGHT - CLIMB TO INITIAL CRUISE ALTI TUDE (LBS) 

WFCR FUEL WEIGHT - CRUISE (LBS) 

WFDES FUEL WEIGHT - DESCEND FROM FINAL CRUISE ALTITUDE (LBS) 

WFRES RESERVE FUEL WEIGHT - TOTAL (LBS) . 

WFRES1 RESERVE FUEL WEIGHT - LOITER + AIR MANEUVER (LBS) 

WFRES2 RESERVE FUEL WEIGHT - ALTERNATE DESTINATION - TOTAL (LBS) 

WFTOT FUEL WEIGHT - TOTAL (LBS) 

WGTO GROSS TAKEOFF WEIGHT (LBS) 

WGULT WGTO X ULTIMATE LOAD FACTOR (LBS) 

WH20 TAKEOFF WATER INJECTION FLUID WEIGHT (LBS) 

WPASS PASSENGER WEIGHT (LBS) 

WPDIST PROPULSION SYS WEIGHT - FUEL SYS DISTRIBUTION (LBS) 

WPDUMP PROPULSION SYS WEIGHT - FUEL SYS DUMPING (LBS) 

WPEC PROPULSION SYS WEIGHT - ENGINE CONTROLS (LBS) 

WPENG PROPULSION SYS WEIGHT - ENGINES (LBS) 

WPFS PROPULSION SYS WEIGHT - FUEL SYS TOTAL (LBS) 

WPL PAYLOAD WEIGHT (LBS) 

WPPUMP PROPULSION SYS WEIGHT - FUEL SYS PUMPS (LBS) 

WPREFL PROPULSION SYS WEIGHT - FUEL SYS REFUELLING (LBS) 

WPROP PROPULSION SYS WEIGHT - TOTAL (LBS) 

WPSEAL PROPULSION SYS WEIGHT - FUEL SYS SEALING (LBS) 

WPSP PROPULSION SYS WEIGHT - NACELLE SOUND PROOFING (LBS) 

WPST PROPULSION SYS WEIGHT - STARTING SYSTEM (LBS) 

WPSYSC PROPULSION SYS WEIGHT - FUEL SYS CONTROLS (LBS) 

WPVENT PROPULSION SYS WEIGHT - FUEL SYS VENTING (LBS) 

WPWI PROPULSION SYS WEIGHT - WATER INJECTION PROVISIONS (LBS) 

WRESCR RESERVE FUEL WEIGHT - ALTERNATE DESTINATION - CRUISE (LBS) 

WRTO WEIGHT RATIO W0/W1 FROM BREGUET EQUATION 

WSB STRUCTURAL WEIGHT - BODY (LBS) 

WSEAC SYS + EQUIP WEIGHT - AIR CONDITIONING (LBS) 

WSEAI SYS + EQUIP WEIGHT - ANTI-ICE (LBS) 

WSEAPU SYS + EQUIP WEIGHT - AUX POWER UNIT (LBS) 

WSEAUX SYS + EQUIP WEIGHT - AUX GEAR (LBS) 

WSEAVI, SYS + EQUIP WEIGHT - AVIONICS INSTALLATION (LBS) 

WSEFRN SYS + EQUIP WEIGHT - FURNISHINGS (LBS) 

WSEHYD SYS + EQUIP WEIGHT - HYDRAULICS ♦ PNEUMATICS (LBS) 

WSEINS SYS + EQUIP WEIGHT - INSTRUMENTS (LBS) 

WSELEC SYS + EQUIP WEIGHT - ELECTRICAL (LBS) 

WSESC SYS + EQUIP WEIGHT - SURFACE CONTROLS (LBS) 

WSH STRUCTURAL WEIGHT - HORIZONTAL TAIL . (LBS) 

WSLG STRUCTURAL WEIGHT - 6ANDING GEAR (LBS) 

WSNAC STRUCTURAL WEIGHT - NACELLES (LBS) 

WSTEW WEIGHT OF STEWARDESS (LBS) 

WSTR STRUCTURAL WEIGHT - TOTAL LBS) 

WSV STRUCTURAL WEIGHT - VERTICAL TAIL (LBS) 

WSW STRUCTURAL WEIGHT - WING (LBS) 

WSYSEQ SYS + EQUIP WEIGHT - TOTAL (LBS) : v 
WT WORKING VALUE - GROSS WEIGHT (LBS) 

WTO GROSS WEIGHT AT START OF CRUISE (LBS) 

WTRTO RATIO TOTAL FUEL WEIGHT/INITIAL CRUISE GROSS WEIGHT 
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TABLE 3-1. DICTIONARY OF VARIABLES (CONT) 

NAME DESCRIPTION 

WTX WORKING VALUE - GROSS WEIGHT (LBS) 

WT1 GROSS WEIGHT AT END OF CRUISE (LBS) 

WUL USEFUL LOAD WEIGHT - TOTAL (LBS) 

WULCR USEFUL LOAD WEIGHT - CREW (LBS) 

WULEQ USEFUL LOAD WEIGHT - EQUIPMENT (L&s) 

WULOIL USEFUL LOAD WEIGHT - ENGINE OIL (LBS) 

WULUF USEFUL LOAD WEIGHT - UNUSABLE FUEL (LBS) 

WUTO FUEL WEIGHT - WARMUP AND TAKEOFF (LBS) 

X LIMIT WING GLOVE THICKNESS (FT) 

XM TABLE - ENGINE DATA - MACH NUMBER 

XMACH WORKING VALUE - MACH NUMBER 

SUBROUTINE CLOES 


ALT 

ALTI 

CL 

DDIST 

DIST 

DTI ME 

DWFL 

GAMAVG 

SFCAVG 

SFC1 * 2 

TGAM1 »2 

THAVG 

THR1 1 2 

TIME 

VCL 

VT1.2 

WFUEL 

WT 

WTX 

XH1 »2 

XMACH 

XMl,2 


WORKING VALUE - ALTITUDE (FT/1000) 

MAXIMUM CLIMB ALTITUDE (FT/1000) 

LIFT COEFFICIENT 

INCREMENTAL CLIMB/DESCENT DISTANCE (FT) 

CUMULATIVE CLIMB/DESCENT DISTANCE (FT) 

INCREMENTAL CLIMB/DESCENT TIME (HR) 

INCREMENTAL CLIMB/DESCENT FUEL WEIGHT (LBS) 

AVERAGE CLIMB GRADIENT (FT) 

AVERAGE SPECIFIC FUEL CONSUMPTION (LB/LB/HR) 

SPECIFIC FUEL CONSUMPTION AT SEGMENT END POINTS (LB/LB/HR) 
CLIMB GRADIENT AT SEGMENT END POINTS (FT) 

AVERAGE TOTAL NET THRUST (LBS) 

TOTAL NET THRUST AT SEGMENT END POINTS (LBS) 

CUMULATIVE CLIMB/DESCENT TIME (HR) 

MAXIMUM CLIMB AIRSPEED (KEAS) 

TRUE AIRSPEED AT SEGMENT END POINTS <KT) 

CUMULATIVE CLIMB/DESCENT FUEL WEIGHT (LBS) 

CURRENT GROSS WEIGHT (LBS) 

CURRENT GROSS WEIGHT ESTIMATE (LBS) 

ALTITUDE AT SEGMENT END POINTS (FT/IOOO) 

WORKING VALUE - MACH NUMBER 
MACH NUMBER AT SEGMENT END POINTS 


subroutine DOVRL 


CD DRAG COEFFICIENT - TOTAL AT XMACH 

CDI DRAG COEFFICIENT - INDUCED 

CDO DRAG COEFFICIENT - PARASITE 

CDOMCR DRAG COEFFICIENT - TOTAL AT CRUISE MACH 

DELCDI COMPRESSIBILITY INCREMENT AT XMACH - CDI 

DELCDO COMPRESSIBILITY INCREMENT AT XMACH - CDO 

DELM MACH NUMBER INCREMENT BELOW CRUISE MACH NUMBER 

DELMDD CRITICAL MACH NUMBER INCREMENT DUE TO CL 

DL RATIO CD/CL AT XMACH 


SUBROUTINE TANGAM 


FN NET THRUST/ENGINE AT XMACH AND ALT (LBS) 

SFC SPECIFIC FUEL CONSUMPTION AT XMACH AND ALT (LB/LB/HR) 

T GAM CLIMB GRADIENT AT XMACH AND ALT (FT) 

THRUST TOTAL NET THRUST AT XMACH AND ALT (LBS) 
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TABLE 3-1 • DICTIONARY OF VARIABLES (CONT) 
NAME DESCRIPTION 

SUBROUTINE TLU1 

I INDEX VALUE FOR X-VECTOR SEARCH 

N DIMENSION OF X AND Y VECTORS 

Q INTERPOLATION FRACTION 

YOUT OUTPUT VALUE OF Y = F C X > 

YVECT ADDRESS OF Y-VECTOR 

XIN INPUT VALUE OF PARAMETER X 

XVECT ADDRESS OF X-VECTOR 

SUBROUTINE TLU2 

J INDEX VALUE FOR Z-VECTOR SEARCH 

NX DIMENSION OF X-VECTOR 

NY INDEX VALUE FOR X-VECTOR SEARCH 

NZ DIMENSION OF Z-VECTOR 

Q INTERPOLATION FRACTION FOR Z-VECTOR 

XIN INPUT VALUE OF PARAMETER X 

XVECT ADDRESS OF X-VECTOR 

YMAT ADDRESS OF Y-MATRIX (DIMENSION NX X NZ ) 

YOUT OUTPUT VALUE OF Y * F(X*Z) 

Yl *Y2 INTERPOLATED VALUES FROM X-VECTORS 

ZIN INPUT VALUE OF PARAMETER Z ' • ■ 

ZVECT ADDRESS OF Z-VECTOR 
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Table 3-2; Dictionary of Variables 


DOC Module (DOCIN) 


AO 

AFMF 

AFSP 

ANHR 

CAVI . 

CCENG 

CnLH 

CFT 

CMLH 

COT 

CREW 

OHM ■ 

ENGSP 

EOT YPE 

EEC 

FR AO 

FRGD 

GO 

GOMVP 

IEQ .. 

IH AX 

MCF 

MLL 

Mil 

PO 

PR 

TURNT- 

YRS 


AIR DELAY - MINUTES 

airframe _maint_enince coeffjlcient 

AIRFRAME SPARES FACTOR 

IDEALIZED ANNUAL UTILIZATION HOURS 

SHIP SEJ ... COST OF. AVIONICS .- $ .. 

ENGINE COST FACTOR 

COST OF MAINTENANCE DIRECT LABOR S/HOUR 

COST OF FUE.L “ $./10jPO LBS. . ... _ 

MANUFACTURING COST PER HOUR $ 

COST OF OIL = S/LB 

CREW . COST PER HCUR.„. _ .... 

DISTANCES (UP TO 20) NAUTICAL MILES 
ENGINE SPARES FACTOR 

FLAG - Q = SUMMARY EQUATIONS, 1 = DETAILED. EQUAT. 
ENGINE UNIT COST 

FUEL RATE - AIR OELAY - 1000 IBS/MINUTE 

FUEL RATE - GROUND DELAY - 1000 LBS /.MINUTE 

GROUND DELAY - MINUTES 

GROUND MANEUVER FIME/FLIGH1 -MINUTES 

FLAG - 1 = RE AO DOCEQ, _0 = OO NOT _READ DOCEQ 

LAST DISTANCE IN DNM - INTEGER TO 20 

MATERIAL ESCALATION COST FACTOR 

MANUFACTURING L A DOR LEARNING EXPONENT 

MANUFACTURING MATERIAL LEARNING EXPONENT 

PRODUCTION QUANTITY - UNITS 

REFERENCE PRODUCTION QUANTITY 

TURNAROUND TIME (HR) 

OEPF.ECIATION LIFETIME - YEARS 


DOC Module; (DOCEQ) 


CEAQ.Q._ 

CEACL 

CEACM 

CEACR 

CEACS 

CEAIC 

CEAIL 

CEAIM 

Cr AIR 

CEAIS 

CEAPUC 

CEAPUL 

CEAPUM 

CEAPUR 


MATERIAL SCALING EXPONENT AIR-CONDITIONING 
LABOR HRS/L8 AT REF W7 - AIR-CONDITIONING 
MATL C0ST/L8 AT REF WT - AIR-CONDITIONING 
REFERENCE WEIGHT FOR AIR-CONDITIONING 
LABOR SCALING EXPONENT - AIR-CONDITIONING 
MATERIAL SCALING EXPONENT - ANTI-ICING 
LAIOR HRS/L9 AT REF WT - ANTI-ICING 

MATL COST/LB AT REF WT - ANTI -ICING 

REFERENCE WEIGHT FOR ANTI-ICING 
t A90R_ SCALING EXPONENT - ANTI-ICING 
MATERIAL SCALING EXPONENT - AUXILIARY POWER UNIT 
LABOR HRS/LB AT REF WT - AUXILLARY POWER UNIT 

MATL C0ST/L9 AT REF WT - AUXILLARY POWER UNIT 

REFERENCE WEIGHT FOR AUXILLARY POWER UNIT 
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Table 3-2. Dictionary of Variables (Cont.) 


C£ AP4JS' LA30R SCALING. EXPONENT - AUXILLARY POWER UNIT 

„ CEAUXC H4TERIAL SCALING EXPONENT - AUXILLARY GEAR 

CEAUXL LABOR HRS /LB AT REF WT - AUXILLARY GEAR 

CEAUXM MATL COST/L 8 AT REF WT - AUXILLARY GEAR 

CEAUXR REFERENCE WEI GHT_ FOR AUXILLARY GEAR 

CEAUXS LABOR SCALING EXPONENT - AUXILLARY GEAR 

CEHYDC MATERIAL SCALING EXPONENT - HYORAULICS/PNEUM ATICS 

CEHYO.L LA80R HRS/LB AT_REF WT - HYDRAULICS/PNEUMATICS 

CEHYOM MATL COST/L 3 AT REF WT - HYORULICS/PNEUMATICS 

CEHYOR REFERENCE WEIGHT FOR HYDRAULICS/PNEUMATICS 

CEHYOS. LABOR SCALING EXPONENT - HYQRAULICS/PNUMATICS 

CEINC MATERIAL SCALING EXPONENT - INSTRUMENTS 

CEINL LA80R HRS/LB AT REF WT - INTRUMENTS 

CEIiiJL MATL COST/LB AT REF. WT .*» INSTRUMENTS 

CEINR REFERENCE WEIGHT FOR INSTRUMENTS 

CEINS LABOR SCALING EXPONENT - INSTRUMENTS 

CESCC MATERIAL SCALING EXPONENT - SURFACE CONTROLS 

CESCL. LABOR HRS/LB AT REF WT - SURFACE CONTROLS 

CESCM MATL C0ST/L9 AT REF WT - SURFACE CONTROLS 

CESCR REFERENCE WEIGHT FOR SURFACE CONTROLS 

CESCS .. LABOR SCALIN£L. EX PON ENT - SURFACE CONTROLS 

CFURNC MATERIAL SCALING EXPONENT - FURNISHINGS 

CFURNL LABOR HRS/LB AT REF WT - FURNISHINGS 

CFURNM . ... MAIL COST/LB AT.. REEL, WT - .FURNISHINGS. 

CFUPNR REFERENCE WEIGHT FOR FURNI9HINGS 

CFURNS LABOR SCALING EXPONENT ~ FURNISHINGS 

CLECC .. . _ MATERIAL ..SCALING —EXPONENT - ELECTRICAL 

CLECL LABOR HRS/LB AT REF WT - ELECTRICAL 

CLECM MATL COST/LB AT REF WT - ELECTRICAL 

CL.ECBL- REFERENCE MEI.GM.I_F.QR_ELEC.TRICAL. 

CLECS LA80R SCALING EXPONENT - ELECTRICAL 

CPECC MATERIAL SCALING 'EXPONENT - ENGINE CONTROLS 

CPE CL.. . L AB.QR. HRS /LB. AT REF WT .--ENGINE CONTROLS 

CPECM MATL COST/LB AT REF WT - ENGINE CONTROLS 

CPFCR REFERENCE WEIGHT FOR ENGINE CONTROLS 

CP ECS LABOR SCALING. EXPONENT ..-.ENGINE CONTROLS. . 

CPFSC MATERIAL SCALING EXPONENT - FUEL SYSTEM 

CPFSL LABOR HRS/LB AT REF WT - FUEL SYSTEM 

CPFSft .MAIL. GOST /LB AT. REE. WT - FUEL SY.3IEM 

CPFSR REFERENCE WEIGHT FOR FUEL SYSTEM 

CPFSS LABOR SCALING EXPONENT - FUEL SYSTEM 

...CPSPC . MATERIAL SC A L ING ^EXPONENT. --- SOUND. PROOFING _ . . 

CPSPL LABOR HRS/LB AWfe# WT - SOONO PROOFING 
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Table 3-2. Dictionary of Variables (Cont.) 


CPSPM 

CP S PR 

CPSPS 

CPSTC 

C P STL 

CP STM 

CPSTR 

CP STS 

CPWIC 

C p WIL 

C°WIM. 

CPWIR 

CPWIS 

CS3C 

CS3L 

CSBM 

CSBR 

CS3S 

CSHC 

CSHL 

CSHM 

CSHR 

CSHS 

CSLGC 

CSLGL 

CSLGM 

CSLGR 

C3LGS 

CSNAC 

CSNAL 

CSNAM 

CSNAR 

CSNAS 

CSVC 

CSVl 

CS VM. _ 

CSVR 

CSVS 

GSWC 

CSWL 

CSWM 

CSMR 

CSWS 


NATL COST /LB AT REF WT - SOUND PROOFING 

REFERENCE HEIGHT FOR SOUND PROOFING : .. - 

LABOR SCALING EXPONENT - SOUND PROOFING 
MATERIAL SCALING EXPONENT - STARTING 
LABOR HRS* LB AT REF WT - SXARTI_NG_ . . ... 

MATL COST/LB AT REF WT - STARTING 
REFERENCE WEIGHT FOR STARTING 
LABOR. SCALING EXPONENT - STARTING .. _ ...... 

MATERIAL SCALING EXPONENT - WATER INJECTION 
LABOR HRS/LB AT REF WT - WATER INJECTION 
MAIL COST/LB AT .REF. WT r WATER .INJECTION 
REFERENCE WEIGHT FOR WATER INJECTION 
LABOR SCALING EXPONENT - WATER INJECTION 
MATERIAL SCALING. EXPONENT - FUSELAGE.. 

LABOR HRS/LB AT REF WT - FUSELAGE 

MATL COST/LB AT REF WT - FUSElAGE 

REFERENCE WEI.GHT...FO.R FUSELAGE 
LABOR SCALING EXPONENT - FUSELAGE 
MATERIAL SCALING EXPONENT - HORIZONTAL TAIL 
LABOR HRS/L9 AT REF WT - HORIZONT AL _T AI L 

MATL COST/LB AT REF WT - HORIZONTAL TAIL 

REFERENCE WEIGHT FOR HORIZONTAL TAIL 
LABOR .SCALING .EXPONENT - HORIZONTAL TAIL ... 
MATERIAL SCALING EXPONENT - LANDING GEAR 
LABOR HRS/L3 AT REF WT - LANDING GEAR 

MAIL COSX/LB AT REF WT - LANDING GEAR . , . 

REFERENCE WEIGHT FOR LANDING GEAR 
LABOR SCALING EXPONENT - LANDING GEAR 

MATERIAL SCALING EXPjDNFNT - NACELLES 

LABOR HRS/LB AT REF WT - NACELLES 

MATL C0ST/L8 AT REF WT - NACELLES 

reference .height fq£. nacelles :_i. .„. 

LABOR SCALING EXPONENT - NACELLES 
MATERIAL SCALING EXPONENT - VERTICAL TAIL 
LABOR HRS/LB AT REF WT - VERTICAL TAIL 

. M4TL. GJ0ST/.L3 AT ..REF WT.. VERTICAL TAIL 
REFFRENCE WEIGHT FOR VERTICAL TAIL 
LABOR SCALING EXPONENT - VERTICAL TAIL 

MATERIAL SCALING EXPONENT. .T.WIN.G. ,.. 

LABOR HRS/LB AT REF WT - WING 

MATL COST/LB AT REF WT - WING 

RE.FERE N G£ . . W£ I G .HJL FOR H.M 

LABOR SCALING EXPONENT - WING 
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Table 3-3. Dictionary of Variables 


ROI Module 
AF 

AN HR 

OELR 

EPS A 

FAR El 

FARES 

FLM 

I RTF 

ISU8 

HXNF 

T 

TAX 

Z( 1) 

Z( 2) 

Z ( 3 ) 

Z ( 4 ) 

7(5) 

Z (6) 


(ROIIN) 

SYS TEH GROWTH FACTOR 3Y CITY PAIR 
ANNUAL HOURS IN SERVICE 

STEPPING FACTOR IN ROI ITERATION-OEC. p RA CT 10 N 
ERROR FRACTION IN ROI ITERATION-DEC. FRACTION 
FIXED CHARGE PER PASS - $ 

CHARGE PER PASS N HI - $ 

MAX. PERMITTED CITY PAIR LOAD FACTOR-DEC FRACTION 

FLAG - 1 = READ ROUTE DATA, 3 = DO NOT READ 

i = PRINT DATA 0Y CITIES , 2_= PRINT SUMMARY ONLY 

MINIMUM FLIGHTS / 0 AY 

DAILY OPERATIONAL PERIOD fHR) 

INCOME T4X RATE .r DEC. FRACTION..; 

SERVICE COST AND LANDING FEE/lOOO LB - DOLLARS 

PASS. SERVICE COST/MILE t DOLLARS 

RFSERV, SALES , PASS HANDLING/PASS - DOLLARS 
STEWARDESS COST/BLOCK HOUR - DOLLARS 

MEAL COST/PASS HOUR - DOLLARS 

GENERAL AND ADMINISTRATIVE IOC FACTOR- PERCENT 
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SECTION 4 


PROGRAM OPERATION AND INSTRUCTIONS 


This section defines the input requirements and physical deck arrangement for each 
of the main program modules operated individually or in sequence . 

Input is divided into data blocks for each of the main program modules and for the 
program control word KEY. The FORTRAN NAMELIST routine requires that each 
input data block be headed by a control card bearing the name assigned to that block. 
Table 4-1 lists the names assigned to each set of input data. In addition to the 
NAMELIST input deck, the ROI module requires a set of city-pair distance and traffic 
data in formatted style (see Section 5). 

A list of input parameters required by each program module, when run independently, 
appears in Table 4-2. Those parameters generated by one model for sequential input 
to the next are indicated; these variables are then omitted from the succeeding input 
deck when running sequentially. Any parameter read into the program supersedes 
the value previously read in or computed. Once a parameter value has been read in, 
it maintains that value unless modified by subsequent input, therefore, it is unnecessary 
to repeat unchanged data for consecutive runs through the same module. 

Figure 4-1 illustrates the deck setup for a sequential run of all three program modules. 
To run modules individually, each input deck must be preceeded by an NKEY data set 
as seen in Figure 4-2. The parameter COM, appearing in the input list of each 
module, provides for a 60 character run title to be reproduced at the head of each 
module output. 

Table 4-1. NAMELIST Input Labels 


PROGRAM ELEMENT 

SUBROUTINE 

INPUT LABEL 

TEKVAL PROGRAM CONTROL 

- 

NKEY 

SIZE 

RNAM 

SIZN 

DOC 

DOIN 

DOCIN, DOCEQ 

ROI 

ROIN 

ROIIN 
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Table 4-2. Input Requirements 


] 

SIZE 

' ALT 

CWEF 

GNFR 

NCREW 

CDF 

CWTIC 

GNFRF 

NENG 

CGLOVE 

DAISLE 

GWAR 

NSTRWC 

CHVOL 

DHCL 

GWS 

PAX 

CLOCH 

DMCL 

GWSWP4 

RANGE 

CLOCV 

DNAC1 

GWTPR 

TRES 

CLOCW 

DRES 

KACS 

VCL 

COM 

DSTP1 

KAF 

WAVI 

CPAXC 

DSTPC 

KAPU 

WCREW 

CRMACH 

DSTW1 

KARULE 

WENG1 

CSWET 

DSTWC 

KEFIN 

WPASS 

CSWPYL 

FLDENS 

KMATL 

WSTEW 

CTWR 

FNRTG 

KPSP 


CWOL 

FNZL 

KWI 


CWFRG 

Tables: 

FTAPER 

NAISLE 

i 


H 

HH 

SFCCL 

DMACH 

HDEL 

XM 

SFCCR 

DCDCV 

MSIGSR 

FNCL 

SFCIDL 

DC DSC 

VDIVM 

FNCR 

SWP 

DCDICV 

DOC (*From SIZE) 

FNIDL 

✓ 1 

FGAM 

DCDISC 

DOCIN: 

AD 

COT 

GD 

PQ 

AFMF 

CREW 

GDMVR 

PR 

AFSP 

DNM 

IEQ 

♦RANGE 

ANHR 

ENGSP 

IMAX 

♦RATING 

CAVI 

EQTYPE 

MCF 

♦TBI 

CCENG 

♦FBI 

MLL 

♦TRS 

CDLH 

♦FBS 

MML 

TURNT 

CFT 

FEC 

♦NENG 

♦WAVI 

CMLH 

FRAD 

♦NSTEW 

♦WGTO 

♦COM 

FRGD 

♦PAX 

YRS 
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Table 4-2. Input Requirements (Cont'd) 


Optional Weights Input (Computed by SIZE): 


For EQTYPE = 0, WAF 




For EQTYPE = 1, 




WSB 

WSLG 

WPFS 


WSEAC 

WSW 

WPSP 

WSESC 


WSEAI 

WSH 

WPST 

WSEINS 


WSEAUX 

WSV 

WPEC 

WSEHYD 


WSEAPU 

WSNAC 

WPWI 

WSELEC 


WSEFRTJ 

DOCEQ (For EQTYPE = 1 Only): 




CSBR 

CSBS 

CSBL 

CSBC 

CSBM 

CSWR 

CSWS 

CSWL 

CSWC 

. CSWM 

CSHR 

CSHS 

CSHL 

CSHC 

CSHM 

CSVR 

csvs 

CSVL 

csvc 

CSVM 

CSNAR 

CSNAS 

CSNAL 

CSNAC 

CSNAM 

CSLGR 

CSLGS 

CSLGL 

CSLGC 

CSLGM 

CPS PR 

CPS PS 

CPSPL 

CPSPC 

CPSPM 

CPSTR 

CPSTS 

CPSTL 

CPSTC 

CPSTM 

CPECR 

CPECS 

CPECL 

CPECC 

SPECM 

CPWIR 

CPWIS 

CPWIL 

CPWIC 

CPWIM 

CPFSR 

CPFSS 

CPFSL 

CPFSC 

CPFSM 

CESCR 

CESCS 

CESCL 

CESCC 

CESCM 

CEINR 

CEINS 

CEINL 

CEINC 

CEINM 

CEHYDR 

CEHYDS 

CEHYDL 

CEHYDC 

CEHYDM 

CLECR 

CLECS 

CLECL 

CLECC 

CLECM 

CEACR 

CEACS 

CEACL 

CEACC 

CEACM 

CEAIR 

CEAIS 

CEAIL 

CEAIC 

CEAIM 

CEAUXR 

CEAUXS 

CEAUXL 

CEAUXC 

CEATJXM 

CEAPUR 

CEAPUS 

CEAPUL 

CEAPUC 

CEAPUM 

CFURNR 

CFURNS 

CFURNL 

CFURNC 

CFURNM 
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Table 4-2. Input Requirements ( Cont’d) 


ROI (*From SIZE, + from DOC) 


AF 

EPSA 

MINF 

ANHR 

FAREI 

*+NSTEW 

*-(COM 

FARES 

*+PAX 

DELR 

FLM 

*+RANGE 

+DNM 

IRTE 

T 


*IMAX 


+DOC 

ISUB 

TAX 


Formatted City-Pair Traffic Data: 

ROUTE TITLE 

NAMEI NAME 2D A P 


*+TBI 

*+TBS 

+TURNT 

*+WGTO 

+XIB 

+YRS 

Z 


JMAX 

PCTM 
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SECTION 5 


SAMPLE RUN 


Following is a reproduction of a sample ease employing all three program modules 
sequentially. The airplane is designed to carry 195 passengers over a design range 
of 5559 km (3000 nm). Cruise speed is to be 0.9 Mach with an initial cruise altitude 
of 11 km (36,000 feet). The design is to employ supercritical aerodynamics, an 
active control system, and composite construction, but no area ruling. The power 
plant is the Pratt & Whitney STF-429 design utilized in the ATT study; sound proofing 
and water injection are included. The configuration calls for three engines in a DC-10 
arrangement. Domestic fuel reserves are indicated. 

DOC input consists of factors typical of domestic operation. A set of detailed cost 
estimation factors are included. Note that the set of distance values DNM are not 
restricted to constant increments but are chosen to reflect the hyperbolic shape of 
the DOC curve. 

Input to the ROI routine lists a system limit load factor of 55 percent and a 48 percent 
tax rate on profit, both values being input. City-pair traffic data are not listed in the 
input, since they are formatted data (not reproduced by the NAMELIST routine); how- 
ever, the information is reproduced in the detailed output listing when called for. 
Figure 5-1 defines the required input card format for traffic data. - 

In the 1972 Convair ATT study, 1985 city pair traffic data were projected! for 156 
domestic and international routes varying in length from 278 km (150 nm) to over 
9265 km (5000 nm). From the 1968 CAB "Handbook of Airline Statistics," the top 100 
city-pairs ranked by number of passengers and number of passenger-kilometers 
yielded data for a total of 128 domestic routes, including mainland-to -Hawaii. Total 
traffic for each route was reduced to a representative fraction for an "average" trunk 
airline by reference to the number of seat-kilometers operated by each trunk author- 
ized to serve that route. Total traffic for an additional 28 U.S. - international routes 
was estimated from data on total seat- kilometers available, average load factors, and 
relative market share of the major airlines serving each route. All data were then 
projected to 1975, 1980, and 1985 by a uniform annual growth rate of 9 percent in 
revenue passenger -kilometers. A complete summary of traffic data is being provided 
to NASA with this model.- - — 

A representative set of 20 domestic routes was selected for the sample problem to 
iUustrate the variation in economic return with distance and traffic level. Traffic 
data are two-way daily passenger totals. Peak hourly rates (see Figure 5-1) are 
simply daily rates divided by the arbitrary 10-hour (input) daily operational period. 
The requirement for peak hourly traffic data is a holdover from previous versions of 
the ROI module and could be deleted with proper program modification. 
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All NAMELIST input data are listed in the sample run as single entries per card, 
followed by a definitive comment. In more common usage, several parameters 
would be listed per card, omitting the definitions, to reduce the bulk of the input 
data decks. 



Figure 5-1. City- Pair Traffic Data Card Format 
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PSNKEY 
KEY =111 , 
SEND 


PJSIZN 


ALT =36. 

,DEF=55H 

COF =-.00270 

» OEF=55H 

CGL OVt =2.4 

,DEF=55H 

CHVOL =0.60 

,OEF=55H 

CLOCH =0.90 

, DEF=55H 

CLOCV =0.95 

,OEF=55H 

CLOCW =0.55 

t DEF = 55H 

COM =60H TEKVAL TEST 

CPA XC =.85 

,0EF=55H 

CRM ACH=-0 • 90 

t OEF=55H 

CSWET =2.02 

»DEF= 55H 

CSW PYL = 0 • 2 

,OEF=55H 

CTWR =0.275 

,DEF=55H 

CVVOL =0.0866 

,DEF=55H 

CWEF =0.77 

,DEF=55H 

CWFRG =0.4 

,0EF=55H 

CWTIC =0.95 

,DEF=55H 

D AIiSLE-21 . 

,OEF=55H 

DHCL =5. 

, OEF= 55 H 

DMCL =0.05 

y DEF=55H 

DNAC1 =7.17 

,DEF=55H 

ORES = 200, 

,DEF=55H 

DSTPC =34. 

f DEF= 55H 

DSTP1 =38. 

,D£F=55H 

DSTWC = 22. 

y DEF=55H 

OSTW1 =25. 

,DEF=55H 

FL0ENS=6.7 

y DEF= 55H 

FNRTG =40000. 

y DEF= 55H 

FNZL =2.5 

,0EF=55H 

FTAPER=4. 33 

,DEF=55H 

GNFR =3.0 

,OEF=55H 

GNFRF =4.0 

y DEF= 55H 

GWAR =9.0 

y OEF= 55H 

GWS =125.6 

y DEF= 55H 


INITIAL CRUISE ALTITUOE(FT/iOOO.) 

SKIN FRICTION COEFFICIENT 

RATIO OF WING ROOT FILLET CHORO/THEC CHORD AT FUS 
HORZ TAIL VOLUME COEFF 
HORZ LOCATION RATIO TO 30DY LENGTH 
VERT LOCATION RATIO TO BODY LENGTH 
WING LOCATION RATIO TO 900Y LENGTH 
CASE - .90 MACH 195 PASS 3000 NM 
PERCENTAGE OF COACH SEATS 
OESIGN CRUISE HACH NO. 

RATIO WING WETTED AREA TO PLANFORM AREA 
PYLON WETTEO AREA FRACTION OF NAC WETTEO AREA 
INITIAL STATIC T/W ESTIMATE 
VERT TAIL VOLUME COEFF 
INDUCED DRAG EFFICIENCY FACTOR 

LATERAL EXTENT OF WING ROOT FILLET/WING SEMI-SPAN 
RATIO WO/WGTO INITIAL ESTIMATE 
AISLE WIOTH (INCHES) 

CLIM8 ALT INCREMENT (FT/100 0.) 

CLIMB MACH INCREMENT 

WING ENGINE NACELLE GIAM<ES = 1) 

RESERVE FUEL DISTANCE (NM) 

COACH CLASS SEAT PITCH(INCHES) 

FIRST CLASS SEAT PITCH (INCHES > 

COACH SEAT WIDTH(INCHES) 

FIRST CLASS SEAT WI DTH(INCHES) 

FUEL DENSITY ( LB/GAL) 

ENGINE STATIC THRUST RATING (ES=1) 

LIMIT LOAD FACTOR 

LENGTH OF TAPERED FUSELAGE SECTIONS (DIAMETERS) 
NACELLE FINENESS RATIO-WING OR FUSELAGE 
NACELLE FINENESS RATIO-FIN 
WING ASPECT RATIO 
WING LOADING 
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GWSWP4=38. ,0EF=55H WING SWEEP ANGLE AT C/4 

GWTPR =0.3374 ,OEF=55M WING TAPER RATIO 

KACS =1 ,0EF=55H ACTIVE CONTROL SYSTEM (Q=WITHOUT, 1=WITH) 

KAF =1 ,0£F=55H AIRFOIL TYPE (3 = CONVENTIONAL 1 = SUPERCRITICAL 

KAPU =1 ,DEF=55H APU ( 0=WITHOUT, 1 = WITH) 

KARULE=0 ,DEF=55H AREA RULE C 0 = STRAIGHT FUSELAGE, 1 = AREA RULED) 

KEFIN =1 ,OEF=55H FIN MOUNTED ENGINE (DC-10 STYLE) ( 0=WITHOUT , 1= WITH) 

KMATL =1 ,DEF=55H CONSTRUCTION MATERIAL ( 0=L I GHT ALLOY , l=COMPOSITE) 

KPSP =1 ,D£F=55H ENGINE NACELLE QUIETING( 0=WIT HOUT, i=WITH) 

KWI =1 ,DEF=55H WATER INJECTION ( 0 = WITHOUT, 1 = WITH) 

NAI SLE= 2 , OEF= 55H NO OF AISLES 

NCREW =3 , DEF= 55H NO IN FLI6HT CREW 

NENG =3 , DEF=55H NO OF ENGINES 

NSROWC=7 *DEF=55H NO OF SEATS PER ROW IN COACH 

PAK =195. ,DEF=55H NO OF PASSENGERS 

RANGE =3000. ,D£F=55H OESIGN RANGE (NM) 

TRES =il. ,DEF=55H RESERVE FUEL TIME (HR) 

VCL =395. ,DEF=55H CLIMB Q-LIMIT (KEAS) 

WAV I =1245. ,0EF=55H AVIONICS SLACK BOX WEIGHTS 

WCREW =195. ,0EF=55H WEIGHT OF CREW MEMBER 

WEMG1 =8326. ,OEF=55H ENG WT(ES = i) INCLUO-TR, NOZZLE , TAIL , SPIN , WATER INJECT 

WPASS =205. ,OEF= 55H PASSENGER WEIGHT 


WSTEW =130 

• 

, OEF= 55H 

WEIGHT 

OF STEWARDESS 




» 

HH = 0. 

» 10 

• , 20. 

,30. , 

36. , 40. 

, 45. , 





XM = . 

3 , • 

4, .5, 

.6, .7 

, • 8 , • 

9, .95, 

.98, 




FNCL(1) = 


24500. , 

23250. , 

22100. , 

20900. , 

19900. , 

19050., 



FNCL (18) = 


21000., 

20000., 

19150. , 

18200. , 

17500. , 

16800. , 

16250 . , 


FNCL (25)= 

o., 

0,, 

o., 

15100. , 

14750. , 

14450. , 

14300., 

14250., 

14200., 

FNCL (37 ) = 

0., 

0.» 

0., 

11520. , 

11510. , 

11510. , 

11500., 

11490 . , 

11490., 

FNCL (49)= 

0., 

0., 

0., 

0. , 

9340. , 

949 0. , 

9640., 

9710., 

9808., 

FNCL (61 ) = 

0., 

0., 

0., 

o., 

7650. , 

7300. , 

7950., 

7950., 

7925., 

FNCL (73)= 

0., 

0., 

0., 

0.1 

5350. , 

6000. , 

6050., 

6100. , 

6100., 

FNCR(l) =2170 0. , 

20400., 

19300., 

18100. , 

17200. , 

163 0 0., 




FNGR(13)=18100. 

, 17300. 

, 16600. 

, 15950., 

15200., 

14600. 

, 14100., 



FNCR(2S)= 

0., 

0., 

o., 

13600. , 

13100. , 

12700. , 

12450. , 

12400., 

12350., 

F NC R ( 3 7 ) = 

0., 

o., 

o.. 

10630. , 

10530. , 

10 480. , 

10350., 

10280. , 

10200., 

F NC R ( 49) = 

0 • , 

o.. 

0., 

0. , 

8780. , 

8750. , 

8850., 

6320 . , 

8800., 

FNCR(6l) = 

0., 

0., 

0., 

o., 

7225. , 

7200. , 

7150., 

7150. , 

7150., 

FNCR(73)= 

0., 

0., 

0., 

0., 

5600. , 

5535. , 

5520., 

5500., 

55.0 0 . , 

SFCCL(l)= 

.560, 

.610, 

. 662 , 

.723, 

.780, 

.884, 
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SFCCL 113) =.540, 

. 585, 

.631, 

.678, 

. 726, 

. 778, 

. 830, 



SFCCLC?5) = 0. , 

0. , 

0. , 

.645, 

.685, 

.729, 

. 770, 

. 791, 

.803, 

SFCCL(37)=0. , 

0. , 

o. , 

.636, 

.671, 

.699, 

.735, 

.750, 

.766, 

S FCCL ( 4 9) = 0. , 

0. , 

0. > 

o. » 

. 666, 

• 694, 

. 720, 

. 731, 

.743, 

SFCCL C 61) = 0. , 

0. , 

0. , 

o. , 

.669, 

. 700, 

.725, 

. 735, 

.753, 

SFCCL C 73) = 0. , 

0. , 

0. , 

o. , 

• 664, 

. 703, 

.730, 

. 747, 

.759, 

SFCCRC 1 ) = .552, 

.615, 

.676, 

.739, 

.803, 

.880 , 




SFCCRC13) = .543, 

.590, 

.639, 

.685, 

.739, 

.795, 

. 850, 



SFCCR(?5) = 0. , 

o. , 

o. , 

. 643, 

.687, 

. 737 , 

. 775, 

.795, 

.610, 

SFCCRC 37) = 0. , 

o. , 

o. , 

.620 , 

.654, 

. 692 , 

. 733, 

.757, 

.765, 

SFCCRC 49) = 0. , 

o. , 

0. , 

o. , 

.655, 

.6 80 , 

.715, 

.732, 

.741, 

SFCCRC 6 1) = 0. , 

0. , 

o. , 

0. i 

. 667 , 

. 700, 

. 721, 

.739, 

.749, 

SFCCRC 73) = 0. , 

0. , 

0. , 

0. , 

.660 , 

.695 , 

. 725, 

. 743, 

.755, 

FNIOLtl)= -7 0 0. 

, -1150., 

-1570. , 

-2120., 

-2670., 

-3220 . , 




FNl DL ( l 3) =-525. 

, —830., 

-1170. , 

-1420., 

-1790., 

-2170 . , 

-2775. 

» 


FNIDL(25)= 0., 

o.. 

o., 

-980. , - 

1225. , 

-1475. , 

-1740., 

-I860., 

-1940., 

FNl DL C 37) = 0., 

0., 

o.. 

-700. , 

-800. , 

-980. , 

-1250. , 

-14 0 35 . , 

-1490., 

FNl OL (49)= fl., 


o., 

o., 

-610. , 

-730. , 

-940., 

-1090. , 

-1190., 

FNl DL C 61 ) = 0-, 

0., 

0., 

0.) 

-500. , 

-60 0. , 

-760., 

-390., 

-980., 

FNII0LI73).- 0., 

0., 

o., 

0., 

-36 0. , 

-450. , 

-600. , 

-700 . , 

-775., 

SFCIPL|1)= 2.0, 

1.28, 

. 879i 

•MV 

• 5 

.506, 



o 

SFC IOLC13)=l. 96 

, 1 *23 , 

.863, 

.732, 

. 604, 

. 521, 

.433 

» 

« . • 

SFCIDL C25) = 0., 

o., 

0., 

. 806, 

.670, 

.5 96 , 

.529, 

. 517, 

.505, 

SFC IDL (37) = 0., 

0., 

O.f 

. 829, 

. 750, 

.654, 

. 552, 

.515, 

.496, 

SFCIDL (49)= 0., 

0.-, 

0., 

o., 

. 304, 

.727, 

. 607, 

. 54 1 , 

.512, 

SFC IDL (61) * 0«, 

0*, 

0 • , 

0., 

. 360 , 

.767, 

.645, 

.585, 

.561, 

SFCIDL (73) - 0., 

0., 

0 . , 

0., 

1.00, 

. 8 66 , 

.684, 

. 615, 

.568, 

VDIVM = 661.7, 

638.1, 

615. 0, 

592.4, 

570.3, 

548.7 , 

527.7, 

507.2, 


487.2, 

467.6, 

448.6, 

430.0, 

412.0, 

394 .4, 

377.2, 



360.6, 

344 .4, 

32 8.6, 

312.7, 

298.7, 

284.7, 

271.4, 



258.6, 

246.4, 

234.9, 







HOEL = 1.0, 

.9298, 

.8637, 

. 3014, 

. 7428, 

.6877, 

.6360, 

.5875, 


.5420, 

.4994, 

.45 95, 

.4223, 

. 3 8 76 , 

.3552, 

.3250, 



.2970, 

.2709, 

.2467, 

.2234, 

.20 38, 

.1351, 

.1681, 



.1527, 

.1387, 

.1260, 







HSI*GSR» ; 1.0, 

.9710, 

.9424, 

. 9143, 

.8 866, 

.8593, 

.3326, 

.8062, 


.7804, 

.7549, 

.7299, 

.7054, 

.6813, 

.6576, 

.6344, 



.6117, 

.5893, 

.5674, 

.5450, 

.5206, 

.4962, 

.4729, 



.4507, 

.4295, 

.40 94, 







H *■ 0., 

2. , 4. 

, 6. , 

8. , 10 

. , 12. 

» 14. , 

16 . , 

18., 20. 

, 22.* 

24. , 

26., 28. 

, 3 0., 

32., 34 

• , 36. 

, 3 8., 

40 . , 

42., 44. 

, 46., 

43. , 











OMACH = 

0. , 
• 45 , 

.01, *02, .04, 

.50, .60, 

.07, 

. 10, 

. 15, 

. 20 

, .25, .30, .35, .40, 

DCOCV = 

0. , 
-11. , 

-4.0, -7. 0, -12., 
-9. 5, -6.0, 

-16. , 

-18., -20., 

-18. 

, -17. ,-15. 5, -14., -12., 

DCDICV= 

0 • 9 

14., 

— 5.0, — 7.0, — 7.5, 
— 1*, . 5 , — 16 . , 

— 8.0, 

—8 .0, —9.0, 

-10. 

,-10. 5, -11. 5, - 12. , -13. , 

DCOSC = 

0. , 
9.5, 

—3.0, —5.0, —9.0, 
— 6.0, 0. , 

-14. , 

— 17., —2 0 . , 

-19.5 

,-18.5,-17.0,-15.0, -12., 

DCD ISC = 

0. , 
-7.0 , 

-1.5, -2.0, -2.0, 
— 8.0, — 10.0, 

-2.0, 

-4.0, -4.5, 

-4.0 

, — 2.0, 0.0, — 3.0, — 5.0, 

SW? 

0. , 

10.0, 15.0, 20.0, 

25.0, 

30 .0, 35.0, 

40.0 

, 45.0, 50.0, 

FSAM = 

.5 , 

.5 , .55 , .62 , 

.74 , 

.89 , 

1.1, 

1.4, 

1.78, 2.26, 


SEND 


en 

i 

05 
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COMMERCIAL TRANSPORT SIZING PROGRAM 

I 

TEKVAL TEST CASE - .90 MACH 195 PASS 3000 NN 
DESIGN SPECIFICATIONS | 


PASSENGERS 

r 

195 

COACH 

165 


FIRST CLASS 

30 


RANGE (NM) 


300 0 

INITIAL CRUISE CONDITIONS 
MACH NUMBER 

l 

. 90 0 

ALTlTUOE( FT) 

1 

36000 

NUMBER OF ENGINES 

i 

3 

MATERIAL 

! 

COMP 

AIRFOIL 


S/C 

ACS 


YES 

AREA RULED FUSELAGE 

i‘ 

1 

NO 


I OERIVEG [SPECIFICATIONS 
GEOMETRY 


FUSELAGE 


LENGTH (FT) 1 

161.63 

DliAMErER(FT) 1 

WING 

17.97 

AREA(S3 FT) i 

1945.23 

T/C i 

. 142 

SPAN (FT) i 

132.31 

ROiOT CHORD (FT) 

21.19 

TIP CHORO(FT) 

8.21 

M A C (FT) 

15.56 

TAIL (EXPOSED AREAS, SQ FT) 


HORIZONTAL 

323.19 

VERTICAL 

178.46 

NACELLES -i 


DIAMETER (FT) i 

5.31 

LENGTH ( FT) ! 

15.93 


WING GEOMETRY 

ASPECT RATIO 9.00 

TAPER RATIO .3874 

SWEEP AT C/4 ( OEG) 38.00 

FUEL RESERVES 

LOITER TIME (HR) 1.00 

OISTANCE TO ALTERNATE (NM) POO 

CRUISE EFFICIENCY AT CRUISE MACH .770 

CDF AT CRUISE MACH . 00270 

.WING LOADING (PSF) 125.6 


PERFORMANCE 

CDOEFF AT CRUISE MACH 

. 01917 

CRUISE CL 

.4368 

CRUISE L/D 

15.637 

ENGINE RATING-LB 

21970 

CLIMB 

DISTANCE (NM) 

511.0 

TIME (MIN) 

60.2 

DESCENT 


DISTANCE (NM) 

75.8 

TIME (MIN) 

11. 1 

BLOCK TIME (HR) 

5.96 

SLOPE (HR/NM) (X1000) 

1.935 

INTERCEPT (HR) (X 1000) 

152.6 

BLOCK FUEL 


SLOPE CLB/NM) 

: 17.864 

INTERCEPT (LB) 

7007.6 



COMMERCIAL TRANSPORT SIZING PROGRAM 


TEKVAL TEST CASE - .90 MACH 195 PASS 3000 NM* 


WEIGHT STATEMENT 


WEIGHT EMPTY 
STRUCTURE 
FUSELAGE 
WING 

HORIZONTAL 
VERT ICAL 
N ACE Ll E S 
LOG GEAR 
PROPULSION 
ENGINES 

SOUND SUPPRESSION 
STARTING 
ENGINE CONTROLS 
FUE. SYSTEM 
PUMPS 

DISTRIBUTION 

'VENTING 

CONTROLS 

REFUEL 

OJMP 

SEALING 

WATER INJECTION 
SYSTEMS AND EQUIPMENT 
SURFACE CONTROLS 
INSTRUMENTS 
HYQR AND PNEU 
AVIONICS 
ELECTRICAL 
AIR CONDITIONING 
AUXILIARY GEAR 
ANTI-ICE 
FURNISHINGS 
ARU 


134.7 
537. 3 

229.4 

119.5 
107.1 

69.2 

312.6 


117199.6 

53372.4 

26 636.1 i 
19271.0 

1740.9 i 

1088.2 
3397.4 ' 

11238.8 

15934.0 

12824.2 
926.9 
16 i. 3 
234.5 

1559.9 


224.3 

■ 37893.2 

3330.7 
819.8 

1224.0 

1618.5 

3101.7 

3463.8 
55.3 

794^3 

22314.5 

1170.6 / 
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7528.1 


USEFUL LOAO 
C RE W 

UNUSABLE FUEL 

ENGINE OIL 

PASS SERVICE EQUIP 

"OPERATING WEIGHT EMPTY 

PAYLOAD 


7 ERO FUEL WEIGHT 

FUEL (TOTAL) 

BLOCK 

WARMUP AND TAKEOFF 
CLIMB 
CRUISE 
DESCENT 
RESERVE 


1365.0 
48 . 6 

164.5 

5950.0 

124727.7 

39975.0 


164899.0 

76853.8 

60600. 2 

1001.3 

16191.7 

43109.9 

296.7 

16253.6 


WATER 


795.9 


GROSS TAKEOFF WEIGHT 


242845.4 
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OI-S 


PSOOCIN 

AO =6. ,DEF=55H AIR OELAY - MINUTES 

AFMF = 1. ,OEF=55H AIRFRAME MAINTENANCE COEFFICIENT 

AFSP = .10 ,OEF=55H AIRFRAME SPARES FACTOR 

ANHR = 4500 • ,OEF=55H IOEALIZEO ANNUAL UTILIZATION HOURS 

CAVI = 500000. ,OEF=55H SHIP SET COST OF AVIONICS - f 

CCENG = 5500. ,OEF=55H ENGINE COST FACTOR 

COL H = 5. ,OEF= 55H COST OF MAINTENANCE DIRECT LABOR $ /HOUR 

CFT = 15. ,OEF=55H COST OF FUEL - $71000 LBS 

.CMLH = 7. ,DEF=55H MANUFACTURING COST PER HOUR $ 

COT- = .66 ,OEF=55H COST OF OIL * f/LB 

CREH = 160. ,DEF=55H CREW COST PER HOUR 

ONM = 0. ,OEF=55H DISTANCES (UP TO 20) NAUTICAL MILES 

ONM = 100., 200., 300., 500., 700., 1000., 1500., 2000., 2500., 3000., 

ENGSP =• .40 ,0EF=55H ENGINE SPARES FACTOR 

EQTYPE= 1 ,0EF=55H FLAG - „ J = SUMMARY EQUATIONS, 1 = DETAILED EOUAT . 

FEC = 0. ,0EF=55H ENGINE UNIT COST 

FRAO =1. ,0EF=55H FUEL RATE - AIR DELAY - 1000 L3S/MINUTE 

FRGO = 1. ,OEF=55H FUEL RATE - GROUND DELAY - 1000 LBS/MINUTE 

GD = 0. , DEF=55H GROUND OELAY - MINUTES 

GDMVR = 0. , OEF= 55H GROUND MANEUVER TIME/FLIGHT-MINUTES 

IEQ = 1 ,OEF=55H FLAG - 1 = READ DOCEQ , 0 = 00 NOT READ OOCEQ 

IMAX =10 , OEF=55H NUMBER OF ONM ENTRIES - MAX OF 20 

MCF =1. ,DEF=55H MATERIAL ESCALATION COST FACTOR 

MLL = .8 ,DEF=55H MANUFACTURING LABOR LEARNING EXPONENT 

MML = .95 ,DEF=55H MANUFACTURING MATERIAL LEARNING EXPONENT 

PQ = 250. ,DEF=55H PRODUCTION QUANTITY - UNITS 

PR = 250. ,0EF=55H REFERENCE . PRODUCTION QUANTITY 

TURNT = .5 ,OEF=55H TURNAROUND TIME - HOURS 

YRS =15. t OEF= 55H DEPRECIATION LIFETIME - YEARS 

SEND 



TT-S 


3 SDOCEQ 
CS3R ; = 

24870. 

,DEF=55H 

CSWR = 

20070. 

,DEF=55H 

CSHR : = 

2450. 

,0EF=55H 

CSVR : = 

2000. 

, DEF=55H 

CSMAR = 

3410. 

,0£F=55H 

CSLGR = 

,1300 0 . 

,DEF=55H 

CPS PR = 

700. 

,DEF=55H 

CPSTR = 

140. 

,0EF=55H 

CPECR = 

195. 

,DEF=55H 

CPWIR = 

260. 

f 0EF=55H 

CPF SR = 

1.6 

,DEF=55H 

CESCR = 

3400. 

,DEF=55H 

CEINR = 

1740. 

,OEF=55H 

CEHYDR= 

2000. 

,DEF=55H 

CLECR = 

3217. 

,0EFs55H 

CEACR = 

3490. 

,OEF=55H 

CEAIR = 

500. 

,DEF=55H 

CEAUXR= 

85. 

,DEF=55H 

CEAPUR= 

90 8. 

,OEF»55H 

CFURNR* 

23169. 

»OEF= 55H 

CS3S = 

.8 

,0£F=55H 

csws , = 

. 8 

t DEF=55H 

CSHS; = 

.8 

,OEF=55H 

csvs i = 

.8 

,DEF=55H 

CSNAS = 

• 9 

,.0EF=55H 

CSLGS = 

1.0 

,DEF=55H 

CPSPS = 

l. 0 

»OEF=55H 

CPSTS = 

1.0 

,0EF=55H 

CPE CS * 

1.0 

,0EF=55H 

CPMIS = 

1. o 

,0EF=55H 

CPFSS = 

.9 

y DEF = 55H 

CESCS = 

.8 

» DEF= 55H 

CEINS * 

1.0 

9 OEF=55H 

CSNYOS= 

. 9 

,0EF=55H 

CLECS = 

. 8 

9 OEF= 55H 

CEACS = 

.8 

,0EF=55H 

CEAIS = 

. 8 

,DEF=55H 


i 


i 

i 


i 

i 

i 


i 

REFERENCE WEIGHT FOR FUSELAGE 

REFERENCE WEIGHT FOR WING 

REFERENCE WEIGHT FOR HORIZONTAL TAIL 

REFERENCE WEIGHT FOR VERTICAL TAIL 

REFERENCE WEIGHT FOP NACELLES 

REFERENCE WEIGHT FOR LANDING GEAR 

REFERENCE WEIGHT FOR SOUND PROOFING 

REFERENCE WEIGHT FOR STARTING 

REFERENCE WEIGHT FOR ENGINE CONTROLS 

REFERENCE WEIGHT FOR WATER INJECTION 

REFERENCE WEIGHT FOR FUEL SYSTEM 

REFERENCE WEIGHT FOR SURFACE CONTROLS 

REFERENCE WEIGHT FOR INSTRUMENTS 

REFERENCE WEIGHT FOR HYDRAULICS/PNEUMATICS 

REFERENCE WEIGHT FOR ELECTRICAL 

REFERENCE WEIGHT FOR AIR-CONDITIONING 

REFERENCE WEIGHT FOR ANTI-ICING 

REFERENCE WEIGHT FOR AUXILLARY GEAR 

REFERENCE WEIGHT FOR AUXILLARY POWER UNIT 

REFERENCE WEIGHT FOR FURNISHINGS 

LABOR SCALING EXPONENT - FUSELAGE 

LABOR SCALING EXPONENT - WING 

LABOR SCALING EXPONENT - HORIZONTAL TAIL 

LABOR SCALING EXPONENT - VERTICAL TAIL 

LABOR SCALING EXPONENT - NACELLES 

LABOR SCALING EXPONENT - LANDING GEAR 

LABOR SCALING EXPONENT - SOUND PROOFING 

LABOR SCALING EXPONENT - STARTING 

LABOR SCALING EXPONENT - ENGINE CONTROLS 

LABOR SCALING EXPONENT - WATER INJECTION 

LABOR SCALING EXPONENT - FUEL SYSTEM 

LABOR SCALING EXPONENT - SURFACE CONTROLS 

LABOR SCALING EXPONENT - INSTRUMENTS 

LABOR SCALIN6 EXPONENT - HY DRAULICS/PNUMATICS 

LABOR SCALING EXPONENT - ELECTRICAL 

LABOR SCALING EXPONENT - AIR-CONOITIONING 

LABOR SCALING EXPONENT - ANTI-ICING 

i 

I . 

i ' 

I 

I 

i 

I 

I- ' 

I 



tn 

i 

to 



LABOR SCALING EXPONENT - AUXILLARY GEAR 

LABOR SCALING EXPONENT - AUXILLARY POWER UNIT 

LABOR SCALING EXPONENT - FURNISHINGS 

LABOR HRS/LB AT REF WT - FUSELAGE 

LABOR HRS/LB AT REF WT - WING 

LABOR HRS/LB AT REF WT - HORIZONTAL TAIL 

LABOR HRS/LB AT REF WT - VERTICAL TAIL 

LABOR HRS/LB AT REF WT - NACELLES 

LABOR HRS/LB AT REF WT - LANOING GEAR 

LABOR HRS/LB AT REF WT - SOUND PROOFING 

LABOR HRS/LB AT REF WT - STARTING 

LABOR HRS/LB AT REF WT - ENGINE CONTROLS 

LABOR HRS/LB AT REF WT - WATER INJECTION 

LABOR HRS/LB AT REF WT - FUEL SYSTEM 

LABOR HRS/LB AT REF WT - SURFACE CONTROLS 

LABOR HRS/LB AT REF WT - INTRUMENTS 

LABOR HRS/LB AT REF WT - HYDRAULICS/PNEUMATICS 

LABOR HRS/LB AT REF WT - ELECTRICAL 

LABOR HRS/LB AT REF WT - AIR-CONDITIONING 

LABOR HRS/LB AT REF WT - ANTI-ICING 

LABOR HRS/L8 AT REF WT - AUXILLARY GEAR 

LABOR HRS/LB AT REF WT - AUXILLARY POWER UNIT 

LABOR HRS/LB AT REF WT - FURNISHINGS 

MATERIAL SCALING EXPONENT - FUSELAGE 

MATERIAL SCALING EXPONENT - WING 

MATERIAL SCALING EXPONENT - HORIZONTAL TAIL 

MATERIAL SCALING EXPONENT - VERTICAL TAIL 

MATERIAL SCALING EXPONENT - NACELLES 

MATERIAL SCALING EXPONENT - LANOING GEAR 

MATERIAL SCALING EXPONENT - SOUND PROOFING 

MATERIAL SCALING EXPONENT - STARTING 

MATERIAL SCALING EXPONENT - ENGINE CONTROLS 

MATERIAL SCALING EXPONENT - WATER INJECTION 

MATERIAL SCALING EXPONENT - FUEL SYSTEM 

MATERIAL SCALING EXPONENT - SURFACE "CONTROLS 

MATERIAL SCALING EXPONENT - INSTRUMENTS 

MATERIAL SCALING EXPONENT - HYORAULICS/PNEOMATICS 

MATERIAL SCALING EXPONENT - ELECTRICAL 
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CEACC = 
CEAIC = 
CEAUXC= 
CEAPUC= 
CFURNC = 
CS9M = 
CSWM = 
CSWM = 
CSVM = 
CSNAM = 
CSLGM = 
CPSPM = 
CP’S T H T = 
CPE CM = 
CPWIM = 
CPF SM = 
CESCM = 
CEINM = 
CENYOMs 
CLECM = 
CEACM = 
CEAIM = 
CEAUXM= 
CEAPUNs 
CFU RNM= 
SEND 


0.9 

,DEF=55H 

MATERIAL SCALING 

1 EXPONENT 

- AIR-CONDITIONING 

0.9 

t DEF= 55H 

MATERIAL SCALING 

i EXPONE 

NT 

- ANTI-ICING 

1. 0 

y DEF= 55 H 

MATERIAL SCALING 

. EXPONENT 

- AUXILlARY GEAR 

1. 0 

y DEF = 55H 

MATERIAL SCALING 

; EXPONENT 

- AUXILLARY POWER Ul 

1.0 

,OEF=55H 

MATERIAL SCALING EXPONENT 

- FURNISHINGS 

19. 

y D£F=55H 

MATL 

COST/L3 

AT 

REF 

WT 

- 

- JSEL AGE 

18. 

y DEF= 55H 

MAf L 

COST/LB 

AT 

REF 

WT 

- 

WING 

20. 

,DEF=55H 

MATL 

C0ST/L3 

AT 

REF 

WT 

- 

HORIZONTAL TAIL 

20. 

,OEF=55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

VERTICAL TAIL 

13. 

y OEF= 55 H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

NACELLES 

19. 

y DEF = 55H 

MATL 

COST /L 8 

AT 

REF 

WT 

- 

LANDING GEAR 

10. 

,DEF=55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

SOUND PROOFING 

76. 

y OEF= 55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

STARTING 

72. 

y DEF= 55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

ENGINE CONTROLS 

75. 

, DEF= 55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

WATER INJECTION 

54. 

y OEF= 55H 

MATL 

COST /LB 

AT 

REF 

WT 

- 

FUEL SYSTEM 

23. 

y D£F= 55 H 

maTl 

COST/LB 

AT 

REF 

WT 

- 

SURFACE CONTROLS 

300. 

y DEF= 55H 

MATL 

C0ST/L8 

AT 

REF 

WT 

- 

INSTRUMENTS 

30. 

,OEF=55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

HYORULICS /PNEUMATICS 

67. 

, OEFs 55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

ELECTRICAL 

42. 

y D£F = 55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

AIR-CONDITIONING 

42. 

y DEF = 55H 

MAiTL 

COST/LB 

AT 

REF 

WT 

- 

ANTI-ICING 

47. 

,DEF=55H 

MATL 

COST/LB 

AT 

REF 

WT 

- 

AUXILLARY GEAR 

84. 

,OEF=55H 

MAIL 

COST /LB 

AT 

REF 

WT 

- 

AUXILLARY POWER UNIT 

19. 

y QEF= 55H 

MAIL 

COST/LB 

AT 

REF 

WT 

- 

FURNISHINGS 


l 

l 

i 


i 

i 

i 

l 



DIRECT OPERATING COST PROGRAM 
TEKVAL TEST CASE - .90 MACH 195 PASS .3000 NM 

total cost per a/c = s 15295462 


AVERAGE AIRFRAME MANUFACTURING COST 



LABOR 

MAT t L 

TOTAL 

FUSELAGE 

405939 

472530 

878469 

MING 

3 06 577 

360527 

667104 

HORIZONTAL TAIL 

68874 

49000 

117874 

VERTICAL TAIL 

53300 

40000 

9330 0 

NACELLES 

1193 0 6 

44322 

163628 

LANDING GEAR 

45500 

245209 

290709 

STRUCTURE TOTAL 

999496 

1211587 

2211084 

SOUND PROOFING 

4900 

7000 

11900 

STARTING 

. 980 

10640 

116 20 

ENGINE CONTROLS 

1365 

14040 

15405 

WAFER INJECTION 

6006 

19500 

25506 

FUEL SYSTEM 

36 

172 

2 08 

PROPULSION RELATED TOTAL 

13287 

51352 

64639 

SURFACE CONTROLS 

177766 

78039 

255805 

INSTRUMENTS 

0 

5220 00 

522000 

Y ORAULI CS/PNE UMATICS 

21327 

1450 35 

166362 

electrical 

73772 

213971 

287743 

AIR CONDITIONING 

17075 

146470 

163545 

ANTI-ICING 

2688 

21995 

24682 

AUXILIARY GEAR 

0 

3995 

3995 

AUXILIARY POWER UNIT 

0 

76272 

76272 

FURNISHINGS 

113528 

440211. 

553739 

SYSTEMS AND EQUIP TOTAL 

406156 

1647988 

2054143 

AIRFRAME HARDWARE TOTAL 

1418939 

2910927 

432 98 66 

ASSEMBLY AND INTEGRATION 

. 


562833 

SUSTAINING TOOLING 



721664 

SUSTAINING ENGINEERING 



379 97 96 

PROFIT AND WARRANTY EXPENSE 



432937 

RECURRING AIRFRAME COST 


• 

12349779 

CRUISE ENGINES (EACH) 



815227 

AVIONICS 



500000 



1SI-S 


OIRECT OPERATING COSTS - ATA EQUATIONS 
AIR DELAY. < HI Nl 6.0 GNO DELAY ( MIN ) 

FUEL RATE, AIR DELAY <1000 LB/MIN) 1.0000 


0.0 


<$ / AIRCRAFT-NM) 


5ND 0 EL AY <10 00 LB/MIN) 


1.0000 


01 ST ANC E 

BLOCK 

BLOCK 

FLIGHT 

flight 

TOTAL 

CREH 

FUEL 

INSURANCE 

MAINTENANCE 

<NH) 

TIKE 

SPEED 

TIME 

SPEED 

- DOC 


+ 01 L 







100.0 

.346 

239.0 .346 

289 

9. 9953 

5957 

2.2644 

.5751 

4.3598 

200. 0 

.540 

370.7 . 540 

3 71 

5.8800 

4644 

1.26 91 

. 3 533 .. 

2.4412 

300.0 

.7 33 

409.3 .733 i 

409 

4. 5031 

4206 

. 93 74 

.2794 

1 . 8017 

500.0 

1. 120 

446.5 1. 120 

446 

3.4106 

3356 

.67 20 

.2202 

1.2900 

700.0 

1.507 

464.6 1.507 

465 

2. 9402 

3706 

.5582 

.1949 

1.0708 

1000.0 

2.087 

479.1 2. 087 , 

479 

2. 5374 

3593 ' 

.47 29 

.1759 

.9063 

1500. 0 

3.055 

491.1 3.055 

491 

2.3131 

3505 

• 40 66 

.1611 

.7784 

2000. 0 

4. 022 

497.3 4.022 j 

497 

2.1759 

3462 

. 3734 

.1537 

.7145 

2500.0 

4. 989 

501.1 4.989 

5 01 

2. 0935 

3435 

.3535 

.1493 

.6761 

3000.0 

5.957 

503.6 5.957 , 

504 

2. 038 7 _ 

341 9 

.3402 

. 1 463 

.6505 

MAINTENANCE COST BREAKDOWN <J / AIRCRAFT-NM) 







DISTANCE 

AIRFRAME 

i 

AIRFRAME ENGINE 

ENGINE MAINT 




,s 


< NM ) 

LABOR 

MATERIAL ' LABOR 

MAT 

ERIAL BURDEN 






100. 0 

• 5064 

l 

.9388 .4651 


7007 1.7487 






20 0.0 

, .2772 

.5077 .2671 


4095 .9797 



t 


1 

300.0 

.2008 

•3640 .2011 


3124 .7234 



DEPRECIATION 

UTILIZ 

50 0. 0 

.1397 

.2490 .1483 


2348 .5183 






700. 0 

.1155 

.1997 .1256 


2015 .4304 






1000.0 

.0938 

.1628 .1087 


1765 .3645 




2. 2003 

1840.6 

1500. 0 

. 07B6 

.1340 .0955 


1571 .3132 




1. 3520 

2335.5 

2000. 0 

. 0709 

.1197 .0889 


1474 . 2876 




1. 0691 

2675.1 

2500. 0 

.0663 

.1111 .0849 


1416 .2722 




. 8428 

3111.0 

3Q00. 0 

. 0633 

.1053 .0823 


1377 .2620 




. 7458 

3378.8 


1 

1 






.6730 

3630.3 



1 






.6164 

3867.0 



1 






.5881 

4002.4 



, i 






.5712 

4090.1 



’ i 






. 5598 

4151.5 


I 

I 


I 



n o n u o u 


PSROIIN 

AF =.33 ,OEF=55H SYSTEM GROWTH FACTOR 3Y CITY PAIR 

DELR =.05 ,OEF=55H STEPPING FACTOR IN ROI ITERATION-DEC. FRACTION 

EPS A, si. 001 , DEF=55H ERROR FRACTION IN ROI ITERATI ON-OEC. FRACTION 

FAREI =•/. ,DEF=55H FIXED CHARGE PER PASS - $ 

FARES =.0644 ,DEF=55H CHARGE PER PASS N MI - 5 

FLM =.55 ,0EF=55H MAX. PERMITTED CITY PAIR LOAD FACTOR-DEC FRACTION 

IRTE = 1, 0EF=55H FLAG - 1 = READ ROUTE DATA, 0 = DO NOT READ 

ISL) 3 =1, DEF=55H 1 = PRINT OATA BY CITIES , 2 = PRINT SUMMARY ONLY 

MINF =4, DEF= 55H MINIMUM FLIGHTS / DAY 

T =10. , DEF=55H OAILY OPERATING PERIOD - MOURS 

TAX =.48 , DEF=55H INCOME TAX RATE - DEC. FRACTION 

Z(l) = SERVICE COST AND LANOING FEE/1000 L9 - OOLLARS 

Z(2> = PASS. SERVICE COST/MILE - DOLLARS 

Z ( 3) = RESERV, SALES , PASS HANOLING/PASS - DOLLARS 

Z C 4) = STEWARDESS COST/BLOCK HOUR - DOLLARS 

ZC5) = MEAL COST/PASS HOUR - OOLLARS 

Z<6> = GENERAL AND ADMINISTRATIVE IOC FACTOR- PERCENT 

Z = 1.37, .0047, .09, 15.03, ,45 , 1.091 

TEND 



RETURN ON INVESTMENT PROGRAM 
TEKVAL TEST CASE - .90 MACH 195 PASS 3000 NM 

CONSTANTS INDEPENDENT OF AIRCRAFT TYPE 

IDEALIZED ANNUAL UTILIZATION, HOURS 4500 

TURNAROUND TIME / TRIP, HOURS .500 

DAILY OPERATIONS TIME, HOURS 10.00 

MAX. PASSENGER LOAD FACTOR .55 

FEDERAL I AX, PERCENT 48.00 

MINIMUM FLIGHTS / DAY / CITY-PAIR 4 

CONSTANTS DEFENDING ON AIRCRAFT TYPE 

PASSENGER CAPACITY 195 

INVESTMENT BASE PER AIRPLANE, $M 17.509 

DESIGN TOGW, L 3 242845 

J. S. OOMESTIC CITY PAIRS - 150 - 250 0 NM 1985 TRAFFIC 

IOC FACTOR, S / 1000 LB TOGW 1.370 

IOC FACTOR, $ / REVENUE PSGR N MI .0047 

IOC FACTOR, $ / PASSENGER .090 

IOC FACTOR, $ / STEWARDESS HR 15.03 

IOC FACTOR, S / TRIP - FOOD .45 

IOC FACTOR, $ t SUBTOTAL IOC $ 1.091 

RANGE, NM 3000 

DEPRECIATION PERIOD, YEARS 15 

PASSENGERS / STEWARDESS RATIO R 

5-17 


CONSTANTS DEPENDING 

ON DISTANCE 




DISTANCE* N MI 

100 

200 

30 0 

500 

BLOCK TIME, WR 

.3460 

.5395 

.7330 

1.1199 

DOC / ACFT-N MI 

9. <3958 

5. 8800 

4. 5081 

3.4106 

TRIP FARE, 1 

13.44 

19.88 

26.32 

39.20 


700 

1000 

1500 

2000 

2500 

3000 

1.5068 

2.0972 

3.0546 

4.0219 

4. 9892 

5. 9565 

2.9402 

2.5874 

2.3131 

2.1759 

2. 0935 

2. 0387 

52.08 

71.40 

103.60 

v;- 

135. 80 

168.00 

200.20 
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TEKVAL TEST CASE - .90 MACH 195 PASS 3000 NM 

U, s. OOMESTIC CITY PAIRS -.150 - 2500 NM 1985 TRAFFIC 


CITY 

PAIR 

DISTANCE 

PAS S/D AY 

HOURLY PEAK 

OOC 

IOC 

NEW YORK 

BALTIMORE 

156 

189 

19 

4799 

2317 

NEW YORK 

BOSTON 

165 

1945 

194 

22950 

12087 

K C 

ST LOUIS 

200 

240 

24 

4704 

2469 

N Y 

BUFFALO 

251 

701 

70 

9102 

4924 

NEW YORK 

NORFOLK 

254 

259 

26 

5221 

2625 

CHICAGO 

MINNEAPOL 

■300 

392 

39 

5410 

2 957 

L A 

S F 

307 

42 0 

42 

5489 

3026 

N Y 

CLEVELAND 

356 

1063 

106 

14955 

7998 

N Y 

CINCINNATI 

504 

376 

38 

6857 

3564 

NEW YORK 

CHICAGO 

626 

1862 

186 

350 91 

18297 

N Y 

ATLANTA 

656 

513 

51 

9983 

5193 

WASHINGTON 

MIAMI 

799 

290 

29 

9025 

4118 

NEW YORK 

N ORLEANS 

1023 

225 

23 

10536 

4366 

NEW YORK 

HOUSTON 

1235 

261 

26 

12145 

50 90 

CHICAGO 

SAN DIEGO 

1502 

13 2 

13 

13 894 

4751 

L A . 

WASH 

1987 

235 

24 

17322 

6708 

S F 

WASH 

2111 

24 9 

25 

18219 

7158 

NEW YORK 

L A 

2129 

10 76 

108 

50459 

24099 

L A 

BOSTON 

2258 

363 

36 

19269 

8 843 

S F 

BOSTON 

2347 

215 

22 

19891 

7319 


INCOME 

LOAD FACT 

PROFIT 

ROI 

fare 

FLT/DAY 

ACFT 

3218 

.2420 

-2 02 7 

-. 1875 

17.05 

4 

.31 

342 74 

.52<+8 

-397 

-. 0103 

17.63 

19 

1.50 

4776 

.3080 

-1246 

-. 1155 

19.88 

4 

,34 

162 36 

.5135 

1150 

. 0605 

23.16 

7 

.65 

6051 

.3321 

-934 

-. 0824 

23.36 

4 

.37 

10314 

.5024 

1013 

. 08 29 

26.32 

4 

.40 

11237 

.53 82 

1416 

. 1103 

26.77 

4 

.40 

318 0 0 

.5449 

4601 

.1295 

29.93 

10 

1.09 

148 44 

.4823 

2300 

. 1329 

39.46 

4 

.53 

88108 

.5305 

18054 

. 1896 

47.31 

18 

2.72 

25287 

.5266 

5258 

. 1923 - 

49.25 

5 

.78 

16976 

.3723 

1993 

. 0907 

58.46 

4 

.71 

16429 

.2890 

794 

. 0 341 

72.88 

4 

.85 

22617 

.3351 

2799 

.0919 

86.53 

4 

.99 

13692 

.1692 

-2575 

-. 0739 

103.73 

4 

1.15 

31746 

.3016 

4013 

. 0894 

134.96 

4 

1.46 

35663 

.3198 

5349 

.10 98 

142.95 

4 

1.54 

155050 

.5016 

41856 

. 2633 

144.11 

11 

4.26 

55327 

.4654 

14152 

. 2175 

152.42 

4 

1.63 

340 27 

.2758 

3545 

. 0702 

158.15 

4 

1.68 
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SYSTEM OAT A 


GROWTH FACTOR 

. 330 

PASSENGERS PER DAY 

11008 

PASSENGER-N MILES/1900 

8550.0 

DOC (i/1000) 

295.3 

IOC (f/1000) 

137.9 

TOC (S/1000) 

433.2 

INCOME (5/1000) 

627.7 

LOAO FACTOR 

.4480 

PROFIT (1/10 00) 

101.1 

RETURN ON INVESTMENT 

.1322 

NUM0ER OF AIRCRAFT 

23.4 
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